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FROM RUBBER BANDS TO RATIONAL MAPS: 

A RESEARCH REPORT 

DYLAN P. THURSTON 

Dedicated to the memory of William P. Thurston, 194-6-2012 

Abstract. This research report outlines work, partially joint with Jeremy Kahn and Kevin 
Pilgrim, which gives parallel theories of elastic graphs and conformal surfaces with boundary. 
On one hand, this lets us tell when one rubber band network is looser than another and, on 
the other hand, tell when one conformal surface embeds in another. 

We apply this to give a new characterization of hyperbolic critically finite rational maps 
among branched self-coverings of the sphere, by a positive criterion: a branched covering 
is equivalent to a hyperbolic rational map if and only if there is an elastic graph with a 
particular “self-embedding” property. This complements the earlier negative criterion of 
W. Thurston. 
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1. Introduction 

This research report is devoted to explaining a circle of ideas, relating: 

• elastic networks (“rubber bands”) and the corresponding Dirichlet energy, 

• extremal length and other “quadratic” norms on the space of curves in a surface, 

• embeddings between Riemann surfaces, conformal or quasi-conformal, and 

• post-critically hnite rational maps. 

The original motivation for this project is the last point, and more specihcally the question 
of when a topological branched self-covering of the sphere is equivalent to a rational map. 
William Thurston hrst answered this question more than 30 years ago [DH93], by giving 
a “negative” characterization: a combinatorial obstruction that exists exactly when the 
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map is not rational (Theorem 8.16). In this paper, we give a “positive” characterization: a 
combinatorial object that exists exactly when the map is rational, for a somewhat restricted 
class of maps (Theorem 1). 

Having both positive and negative combinatorial characterizations for the same property 
automatically gives an algorithm for testing the property. You search, in parallel, for either 
an obstruction or a certihcate for the property. One of the searches will eventually succeed 
and answer the question. There are other algorithms for testing whether a branched self-cover 
is rational: you may take the rational map itself as a certihcate for being rational [BBY12]. 
However, we expect our combinatorial certihcate to be more practical to search for (Section 6.6). 
By contrast, W. Thurston’s obstruction theorem is notoriously hard to apply. 

In addition, a positive characterization gives an object to study associated to the topological 
branched self-covers of most interest, namely the rational ones. 

Several of the constructions along the way are of independent interest. For instance, we 
give a new characterization of when one Riemann surface conformally embeds in another in a 
given homotopy class, and a numerical invariant of such embeddings. 

This is a preliminary report on the results. Most proofs are omitted, although we give some 
indications. We also indicate some of the many open problems suggested by this research 
in the subsections titled “Extensions” (Sects. 3.4, 5.5, 6.7, 7.4, and 8.7); all of these are 
unnecessary for the main results. 

1.1. Detecting rational maps. Recall that a (topological) branched self-cover of the sphere 
is a hnite set of points P in the sphere and a map / : (F^, P) (5^, P) that is a covering when 
restricted to S\f~^{P). One central question is when such a map is equivalent to a rational 
map on CP^. (See Definition 8.1.) A branched self-cover can be characterized (up to homotopy) 
by picking a spine T for S\P, and drawing its inverse image F = /“^(F) c S\f~^{P). 
There are two natural homotopy classes of maps from F to F. (In this paper, a map of graphs 
is a topological map, not necessarily taking vertices to vertices.) 

• A covering map commuting with the action of /: 


S\f-\P) S\P. 

We denote this covering map tt : F —> F. 

• A map commuting up to homotopy with the inclusion in S\P: 


S\f-\P) S\P. 

We denote this map cf. It is unique up to homotopy since F is a spine for S\P. 

Example 1.1. Figure 1 shows a simple example of such spines when P has 3 points. A, R, 
and C (at oo). The covering map tt preserves the colors and arrows on the edges. The 
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Figure 1. A spine and its inverse image defining a branched self-cover of 
the sphere. The marked set P consists of the two points shown, and a third 
point C at infinity. The nnmbers on the edges give a measnre which proves 
this branched self-cover is eqnivalent to a rational map. 


branched self-cover / : (S*^, P) O is an extension of vr to the whole sphere. This extension is 
nnique np to homotopy relative to {A, B, C}. The map / permntes the marked points by 


( 2 ) 

A'^ * B 



(For instance, A is inside a crimson-goldenrod region in F, so must map to B, which is inside 
a crimson-goldenrod region in F.) 

The map 0, on the other hand, is the projection of F onto F considered as a spine; for 
instance, it might map the right-hand purple edge of F to the right-hand goldenrod edge of 
r, whereas vr preserves the colors. 

For our characterization of rational maps, we also need an elastic structure on F, by which 
we mean a measure a on F absolutely continuous with respect to Lebesgue measure.^ We can 
pull back a by TT to get a measure 7r*a = S on F. For : F ^ F a Lipschitz map of graphs, 
define the embedding energy by 

(1.1) Emb(V’) := esssup V 

yGF 

^ xgT 

^{x)=y 

The derivative is taken with respect to the two measures a and a. The essential supremum 
ess sup ignores points of measure zero. In particular, we may ignore vertices of F or points of 
F that map to vertices of F. 

In practice the embedding energy is optimized when is piecewise linear, and the reader 
may restrict to that case. 

For one motivation for Equation (1.1), see Proposition 6.23: the embedding energy charac¬ 
terizes conformal embeddings of thickened versions of the graph. Another one is in Theorem 6: 
it characterizes when Dirichlet energy is reduced under the map, i.e., when one rubber band 
network is looser than another. (See Remark 6.19.) 

Example 1.2. To return to Example 1.1, consider the measure a on F from Figure 1 and the 
concrete map : F ^ F that maps 

• the right purple edge of F to the goldenrod edge on the right of F, 

• the right goldenrod edge of F to the crimson edge in the middle of F, and 

^We can also interpret a as a metric. But we prefer to distinguish this measure from other metrics that 
come later. 
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• the remaining four edges (two crimson, one goldenrod, and one purple) of F map to 
the purple edge on the right of F, with the crimson edges mapping to segments of 
length l/V^. 

Make 'ip linear on each segment described above. Then a short computation shows that 
Emb(V’) = l/\/2. 

A major theorem of this paper is that rational maps can be characterized by maps ip with 
embedding energy less than 1. We say that two branched self-covers /i : (S'^,Pi) O and 
/2 : (5^, P 2 ) O are equivalent if they can be connected by a homotopy of /i relative to Pi 
and conjugacy taking Pi to P 2 . 

Theorem 1. Let f : {S‘^,P) LL> be a branched self-cover of the sphere. Suppose that there 
is a branch point in each cycle in P. Then f is equivalent to a rational map iff there is an 
elastic graph spine (F, a) for S\P, an integer u > 0, and map p! e [0„] so that Emb('?/') < 1- 

Here (pn is obtained from (p by iteration, and is in the homotopy class of the projection 
y-n(F) ^ p Loosely speaking. Theorem 1 says that the self-cover is rational iff there is a 
self-embedded spine for S'^\P. 

Remark 1.3. It is likely the condition on Theorem 1 can be relaxed to assume merely that / 
has at least one branch point in one cycle in P, i.e., if / is rational, its Julia set is not the 
whole sphere. See Section 8.7. 

Example 1.4. The explicit measure and map in Example 1.2 show that the given branched 
self-cover is equivalent to a rational map. In fact, it is equivalent to the rational map 


with A = 0, P = 1, and C = 00 .^ 

There was a previous characterization of rational maps by W. Thurston, recalled as 
Theorem 8.16 below. This is analogous to the two ways to characterize pseudo-Anosov surface 
automorphisms, which form a natural class of geometric elements of the mapping class group 
of a surface. Geometrically, pseudo-Anosov diffeomorphisms are those whose mapping torus 
is a hyperbolic 3-manifold. Combinatorially, there are two criteria: 

• Negative: f is pseudo-Anosov iff it is not periodic (/°^ is not the identity for any 
fc > 0) or reducible (there is no invariant system of multi-curves for /). 

• Positive: f is pseudo-Anosov iff there is a measured train track T and a splitting 
sequence from T to a train-track T' with T' = A/(T) for some constant 0 < A < 1 
[PP87]. 

The positive criterion gives some extra information: the number A is an invariant of /, with 
dynamical interpretations. (For instance, A controls the growth rate of intersection numbers.) 

Analogously, we can say that a branched self-cover / of {S^, P) is geometric if it is equivalent 
to a rational map, which also have associated 3-dimensional hyperbolic laminations [LM97]. 
There are combinatorial criteria for / to be rational: 

• Negative: f is rational iff there is no obstruction, as in W. Thurston’s Theorem 8.16. 
In loose terms, the obstruction is a back-expanding annular system: a collection of 
annuli that get “wider” under backwards iteration. 

^Every branched self-cover with only 3 post-critical points is equivalent to a rational map, so this example 
was easy to do by other means. 
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• Positive: Under some additional assumptions, / is rational iff there is a metric spine 
for S‘^\P satisfying a back-contracting condition. This is Theorem 1. 

As in the case of surface automorphisms, the two theorems are in a sense “dual” to each other: 
it is easy to see that a branched-self cover cannot simultaneously have a back-expanding 
annular system and a back-contracting spine. (See Equation (8.2).) Also as in the surface 
automorphism case, the positive criterion gives us a new obiect to study, namely the constant 
SF[0] of Section 7.2. 

Compared to the situation for surface automorphisms. Theorem 1 has the following caveats: 

• It only works when there is a branch point in each cycle in P. (But see Remark 1.3.) 

• We may need to pass to an iterate to get Emb(0) < 1. Furthermore, it is easy to see 
in examples that the embedding energy can decrease in powers. (In the notation of 
Section 7.2, in general SF[0„] < SF[0]”' and SF[0] < SF[0]. See Example 2.7.) 

There is a further caveat for both surface automorphisms and branched self-covers: 

• For the positive criterion, the train track or graph constructed is not canonical: there 
are many different choices that work for the criterion. 

By contrast, the negative criteria can be made canonical. (Pilgrim [PilOl] proved this for 
branched self-covers.) On the other hand, in the surface case, Agol [Agoll] and Hamenstadt 
[Ham09] give a canonical object related to the measured train track. 


1.2. Elastic graphs. The embedding energy of Equation (1.1) looks a little mysterious; it 
looks a little like the Lipschitz stretch factor, but the sum over inverse images looks unusual. 
To explain where it comes from, we now turn to a “conformal” theory of graphs parallel to 
the conformal theory of Riemann surfaces. The central object is an elastic graph (F, a), which 
you should think of as a network of rubber bands; formally, it is a graph with a measure on 
each edge, representing the elasticity of the edge. (See Section 3.3.) 

There are additional structures we can put on the graph. 

• On one hand, we can consider curves C on the graph, maps from a 1-manifold into F. 

• On the other hand, we can consider maps from F to a length graph K, a graph with 
hxed lengths of edges (like a network of pipes). 

Maps between these objects have naturally associated energies, as summarized in the following 
diagram. 





Curve EL Elastic 

Dir Length 

C graph F 

graph K 

u 

u 

Emb 

Lip2 


The labels on the arrows indicate the type of energy on a map of this type, as follows. 

• For a map / from an elastic graph (F, a) to a length graph {K, i), there is the Dirichlet 
or rubber-band energy (Section 4) familiar from physics: 

Dir(/) := f 

JxeT 


(1.3) 
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where /' measures the derivative with respect to the natural metrics. If / minimizes 
this energy within its homotopy class, it is said to be harmonic. 

• For a curve C in an elastic graph (F, a), we have a version of extremal length (Section 5): 

(1.4) EL[C]:= 21 nc(eY-a(e), 

eGEdge(r) 

where nc{e) is the number of times C runs over the edge e (without backtracking). 

• For a curve C in a length graph {K, i), we have the usual length, which in our notation 

is 

(1.5) ^[C]-= nc{e)-i{e). 

eeEdge{K) 

To match the other quantities, we actually use the square of the length as our energy. 

• For a map 0 between length graphs and (7^2,£ 2 ), there is the Lipschitz 

constant (Section 4.3): 

(1.6) Lip(0) := esssup|0'(a;)|. 

xgKi 

Again, we consider the square of the Lipschitz energy. 

• Finally, for a map 0 between elastic graphs (Fi,a;i) and (F 2 ,q; 2 ), we have embedding 
energy as used in Theorem 1 (Section 6): 

(1.7) Emb(0) := ess sup L m^)i 

This energy appears to be new, although it is related to Jeremy Kahn’s notion of 
domination of weighted arc diagrams (Section 8.5). 

Remark 1.5. We could make the diagram more symmetric by using width graphs instead of 
curves (see Section 3.3), and adding a norm on maps between width graphs. 

These energies are sub-multiplicative, in the sense that composing two maps can only 
decrease the product of the energies: if / and g are two composable maps of the above types, 
then 

(1-8) ll/ll-ll^ll, 

where ||'|| is the appropriate energy from the above list. (This inequality is the reason we 
squared some of the energies.) For instance, if we hx elastic graphs Fi and F 2 , a length 
graph K, and maps 0 : Fi ^ F 2 and / : F 2 —> K, then 

(1.9) Dh(/ o 0) ^ Dh(/) • Emb(0) 

What is more, these inequalities are all tight, in the sense that if we hx the domain, range, 
and homotopy type of /, then we can hnd a sequence of functions Qi (including a choice of 
domain) that approach equality in Equation (1.8). Likewise if we hx g and vary /, we can 
hnd a sequence of functions fi approaching equality in Equation (1.8). 

For instance, if0:Fi^F2isa map of elastic graphs, we can strengthen Equation (1.9) to 

Dir[/ o 0] 
sup TV- rri ’ 

K length graph J 


( 1 . 10 ) 


Emb[0] 
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where Emb[0] and Dir[/] are the minimums over the respective homotopy classes (Theorem 6). 
Since Dirichlet energy can be interpreted as the elastic energy of a stretched rubber band 
network, Emb[0] < 1 can therefore also be interpreted as saying that Ti is “looser” than r 2 , 
however the two rubber band networks are stretched. See Remark 6.19. 

Other examples are in Propositions 5.13 and 5.14. This gives a kind of duality between 
curves in an elastic graph T and maps from T to length graphs. If we think of curves as living 
in a “vector space” and maps to length graphs in its “dual”, then the embedding energy can 
be interpreted as an “operator norm”. 

This theory of conformal graphs is largely parallel to the theory of conformal (Riemann) 
surfaces with boundary, where we again have a number of energies: 





Curve EL Riemann Dir Length 
C surface S graph K 



Again, each arrow is marked by the appropriate energy for measuring that type of map. 
Dir and EL are again Dirichlet energy and extremal length, but on surfaces rather than 
graphs. SF is new; it is the stretch factor of a homotopy class of a topological embedding 
[0] ; El ^ E 2 between Riemann surfaces. In general, we do not know a direct expression 
analogous to Equation (1.1), so SF is dehned to be the minimal ratio of extremal lengths 
(Dehnition 6.2, analogous to Equation (1.10)). When there is no conformal embedding of Ei 
in E 2 in the given homotopy class, there is a direct expression: SF[0] is given by the minimal 
quasi-conformal constant in the homotopy class (Theorem 3). There is an analogue of SF for 
maps between elastic graphs, and in that context SF[0] = Emb[0] (Theorem 6). 

Remark 1.6. To prove Theorem 1, we actually do not need to consider length graphs or the 
Dirichlet energy at all. (The proofs go through extremal length instead.) However, they 
illuminate the overall structure. In particular, it is not clear why one would consider elastic 
graphs without the rubber-band motivation. 

Remark 1.7. The appearance of length squared in (1.2) and (1.11) is easy to justify on the 
grounds of units. Extremal length itself behaves like the square of a length, in the sense that 
if we take k parallel copies of a curve, the extremal length multiplies by k'^. Likewise, if the 
lengths on the target of a harmonic map are multiplied by k, the harmonic map remains 
harmonic while the Dirichlet energy is multiplied by k"^. 

1.3. History and prior work. Although Equation (1.1) appears to be new, Jeremy Kahn’s 
notion of domination of weighted arc diagrams [Kah06] is essentially equivalent. See Sec¬ 
tion 8.5. 

There has also been substantial work in the related setting of resistor networks rather than 
spring networks [BSST40,Duf62,CIM98, inter alia]. See Remark 4.4. 

Theorem 1 is closely related to Theorem 8.4, which characterizes when a rational map 
exists in terms of conformal embeddings of surfaces. Theorem 8.4 has been a folk theorem in 
the community for some time. 
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For polynomials, Theorem 1 reduces to a previously-known characterization in terms of 
expansion on the Hubbard tree; see Theorem 8.11 in Section 8.3. 

1.4. Organization. After a section giving some examples of how to apply Theorem 1, this 
paper is organized by topics moving up a dynamical hierarchy. For each topic we give hrst 
the conformal surface notions and then the graph notions. 

• We start with notions depending only on a single conformal surface or graph: Dirichlet 
(or rubber band) energy of (harmonic) maps (Section 4) and extremal length of curves 
(Section 5). 

• Next come notions depending on a map between surfaces or graphs. This is the 
stretch factor or embedding energy, which generalizes the Teichmuller distance and 
characterizes conformal embeddings of surfaces (Section 6). 

• Next is the dynamical theory of iterated maps (Section 7). Here we hnd another 
number, the asymptotic stretch factor, which characterizes rational maps (Section 8). 

The paper is organized by a logical hierarchy, rather than what is necessarily pedagogically 
best; the reader is encouraged to skip around. 

1.5. Acknowledgements. I would like to thank Matt Bainbridge, Steven Gortler, Richard 
Kenyon, Sarah Koch, Tan Lei, Dan Margalit, and Giulio Tiozzo for many helpful conversations. 
I would like to especially thank Maxime Fortier Bourque, who pointed me towards Ioffe’s 
theorem [Iof75] and had numerous other insights. 

This project grew out of extensive conversations with Kevin Pilgrim, who helped shape 
my understanding of the subject in many ways. Many of the arguments were developed 
jointly with him. Notably, he communicated Theorem 8.4 to me, and Theorem 4 is joint work 
with him. Theorem 5 is joint work with Jeremy Kahn, who also contributed substantially 
throughout. 

Above all, I would like to thank William Thurston for introducing me to the subject and 
insisting on understanding deeply. 


2. Examples 

Here, we give some more substantial examples of Theorem 1. 

2.1. Polynomials: The rabbit and the basilica. Theorem f is not very interesting for 
polynomials, as every topological polynomial with a branch-point in each cycle is equivalent 
to a polynomial. The extension of Theorem 1 to the general topological polynomial case is 
somewhat more interesting, but is equivalent to known results on expansion on the Hubbard 
tree; see Section 8.3. Nevertheless, we will look at some examples, both to see what the 
stretch factors are and to use them for matings. 

Example 2.1. We hrst look at the “rabbit” polynomial, the post-critically hnite polynomial 
fi{z) = + c with c sa —0.1226 -I- 0.7449L The critical point moves in a 3-cycle 
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The optimal elastic graph Ti and its cover Ti are 




22/3 




with Emb[ri ^ Ti] =2 which is less than one, as expected. (There is another marked 
point at inhnity, not shown.) 


Example 2.2. Another graph r 2 that works to prove that the rabbit polynomial is realizable is 


( 2 . 1 ) 




The black edges have the indicated lengths, which come from looking at the external rays 
landing at the a hxed point of /i. Give the colored edges an equal and long elastic length 
(say, 100). There is a natural map 0o : r 2 —»■ r 2 as follows. 

• The outside circle is mapped to the outside circle, with derivative |0o| = 1/2. 

• The colored segments on the lower right of r 2 is squashed out to the lower-right 
boundary of r 2 , with derivative |0 q| on order of 1/100. Thus Filing % 1/2 on the 
corresponding portion of r 2 . (Fill^ is dehned in Equation (4.4), and is the quantity 
maximized in the definition of Emb((/)).) 

• The colored segments in the upper left of r 2 are mapped to the colored segments in 
r2, with |(/g| = 1. 

Because of the last point, this map has Emb[0o] = 1? which is not good enough. Let (f)i 
be the result of “pulling in” very slightly the image of the ends of the upper-left colored 
segments of r 2 , so they map to the interior of the colored segments of r 2 . This decreases the 
derivative on the colored segments to less than one, while only increasing the derivative on 
the outside circle slightly; thus Emb[0i] is very slightly less than 1, as desired. 


Example 2.3. We can perform a similar trick for other polynomials. For instance, the basilica 
polynomial — 1 has a spine Ts, with cover given by 



The same argument as above (giving the purple edge a long length, and pushing the right 
purple edge out to the boundary) shows that there is a map </> ; Ts —> Ta with Emb[0] < f. 
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(In this case the optimal stretch factor is 2 and is not realized with a graph with this 
topology). 

2.2. Matings. We can nse the techniqnes of Section 2.1 to show that some matings of 
polynomials are geometrically realizable. 


Example 2.4. We can glue together the hgures in Equations (2.1) and (2.2): 


(2.3) 




This gives a graph spine r 4 and cover r 4 representing the formal mating of the rabbit and 
the basilica. We can End a map 0 : r 4 ^ r 4 with embedding energy less than 1 : 

• Assign the black mating circle in r 4 total length 1, divided according to the angles of 
the external rays. 

• Give all colored edges an equal and large length. 

• Pull this metric back to r 4 , and map r 4 to r 4 by pushing colored arcs out to the 
black circle as in Examples 2.2 and 2.3. 

• Pull the map in slightly where colored vertices meet the black circle. 

The result has Fill 0 (|/) slightly larger than 1/2 when y is on the black circle, and Fill^(j/) 
slightly less than one on the colored edges. 


Naturally, the technique of Example 2.4 cannot always work, as sometimes the mating is 
not geometrically realizable. 


Example 2.5. If we try to mate a basilica with a basilica, we get these graphs: 


(2.4) 




If we try to use the same technique as before, it doesn’t work, as there are two points on the 
black circle that we attempt to pull in two different directions. Indeed, the left green-purple 
circle is mapped to a green-purple circle, so must have derivative at least 1 : it is an obstruction 
to the mating. (In this case, it is a Levy cycle). 


2.3. Slit maps. Given a branched self-cover / : (A^, P) O and an arc A with endpoints 
in P, there is a blowing up construction which produces a map that agrees with / outside 
of a neighborhood of A and maps that neighborhood surjectively on to 3“^. Pilgrim and Tan 
Lei showed that, if / is a rational map, these blow ups frequently are as well [PT98]. We will 
restrict attention to cases where the initial map / is the identity, in which case the theorem 
becomes the following. 
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Theorem 2.6. Let P ^ be a finite graph, and let G ^ S'^ be a finite embedded graph with 


endpoints on P. Then idc is a rational map iff G is connected. 

We can give a new proof in the harder direction, when G is connected. We start with a 
simple example. From the connected planar graph 



take as spine r 4 the spherical dnal to G. Then r 4 is obtained by taking the connect snm 
of r 4 with fonr extra copies F^ of F 4 , one for each edge: 



For any metric on F 4 , there is a natural map (/> : F 4 ^ F 4 that maps most of each copy F^ to 
a point in the center of the corresponding edge of Fq. This map has derivative equal to 1/2 
or 0 everywhere. We have Emb(0) = 1/2, and consideration of the red edge shows that this 
is optimal. 

This example generalizes immediately to a general connected graph G, except that Emb[0] 
will only be 1/2 when G has a univalent vertex; otherwise, Emb[0] will be strictly smaller. 
We have thus proved Theorem 2.6, with some additional information about the stretch factor. 

2.4. Behavior under iteration. We can iterate Example 1.1: 



(See Section 7.1 for details on iteration.) In this example, the embedding energy behaves well, 
in the sense that the fc’th iterate fik is optimal for embedding energy (using the edge lengths 
and concrete initial map from Example 1.2): 


Emb[0fc] = Emb((/fc) = (V2) ^ 


Such good behavior is not generally the case. 
Example 2.7. The map 
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SF 



1^I_I_I_I^I_L 

012345678 


n 


Figure 2. Logarithmic plot of SF[(;/)„] for 0 from Example 2.7, as n varies. 
The parameters {a,/3,'y) are the elastic lengths of the left, right, and middle 
edges in Eqnation (2.5), respectively. When 7 = 0, the graph Es degenerates 
to a rose graph. The line for (a,/ 3 , 7 ) = (1,1, —1) corresponds instead to a 
barbell graph homotopy eqnivalent to Es. The dashed line shows SF[(/)]"^. 


is represented combinatorially by a graph E 5 and maps tt, 0 : E 5 ^ E 5 : 



A case analysis shows that for any metric on E and any map "0 e [0] with Emb('0) < 1, there 
is some power 1 < /c < 4 so that the fc-fold iterate i/’fc has local back-tracking. This implies 
that 

Emb['^4] < Emb['^]^, 

regardless of the initial metric. Similar facts hold for any of the other three graphs homotopy 
equivalent to E 5 . Thus, regardless of the choice of initial spine, we have SF[(;/)] < SF[(;/)] = 
Emb[0]. (See Section 7.2 for dehnitions.) 

Figure 2 gives a sample of experimental data for this map. An ad hoc argument shows 
that for this map, SF[0] = 2“^/^. The corresponding line is shown dashed in the hgure. 

3. Setting 

3.1. Surfaces. We work with compact, oriented surfaces S with boundary. It is sometimes 
convenient to think about the double DYl = S ug S, which has no boundary. 

Definition 3.1. A curve on E is an immersion of a 1-manifold with boundary in E, with 
boundary mapped to the boundary. The 1-manifold need not be connected; if it is, the curve 
is said to be connected. Curves are considered up to homotopy within the space of all maps 
taking the boundary to the boundary (not necessarily immersions). A curve is simple if it is 
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embedded (has no crossings). An arc is an interval component of a cnrve, and a loop is a 
circle component. A curve of type + is a curve with only loop components and a curve of 
type — is a curve with no loops parallel to the boundary. The geometric intersection number 
i([Ci], [ 6 * 2 ]) of two curves is the minimal number of intersections (without signs) between 
representatives of the homotopy classes [Ci] and [C 2 ]. 

A weighted curve is a positive linear combination of curves, where two parallel components 
may be merged and their weights added. The space of weighted simple curves on S is denoted 
C(S). If S has boundary, then we distinguish two subsets: 

• C+(S) is the space of weighted simple curves of type +; and 

• C^(S) is the space of weighted simple curves of type —. 

Remark 3.2. As curves need not be connected, they are what other authors would call a 
multi-curve. 

There are several different geometric structures one can put on a surface. First, we can 
consider a conformal or complex structure cc on S, considered up to isotopy. 

The next two structures deal with measured foliations (or equivalently measured lamina¬ 
tions). We always consider measured foliations up to homotopy and Whitehead equivalence. 
Given a measured foliation F and a curve G, we can compute f([G], F), the minimal (trans¬ 
verse) length of any curve isotopic to C with respect to F. This is unchanged under Whitehead 
equivalence, and the converse is true; two measured foliations are Whitehead equivalent iff 
the transverse lengths of all curves within an appropriate dual class (as specified below) are 
the same. 

For a surface S with no boundary, there is only one type of measured foliation, and they 
form a hnite-dimensional space fAFiTi). The space of curves C(S) is dense in AfJ^(S), and 
measuring lengths on C(S) gives an embedding i : Af-F(S) 

On a surface with boundary, measured foliations come in two natural flavors. is 

the space of foliations that are parallel to the boundary (i.e., so the transverse length of the 
boundary is 0 ). C+(S) is dense in AfJ^^(S), and measuring lengths of curves in C^(S) gives 
an embedding i : 

Dually, fAF^ifS) is the space of measured foliations without boundary annuli. It is the 
closure of C~(S), and there is an embedding i ; JAF^ifF) ^ 

For a closed surface S, we dehne JAF^iTi) and JAF^ {F) to be equal to JAFfF). 

Warning 3.3. The set of connected curves is not dense in JAF-{F) if E has non-empty 
boundary. 

Finally, we consider quadratic differentials, in the following variants. A quadratic differential 
on S is locally of the form q{z) [dzY with q holomorphic, and determines a half-turn surface 
structure on S, away from a hnite number of singular cone points. (A half-turn surface is a 
surface with a chart where the overlap maps are translations or rotations by tt.) Q(E,a;) is 
the space of quadratic differentials with hnite area. 

If S has boundary, then Q®(S,a;) is the space of quadratic differentials that are real on 
the boundary; this is isomorphic to the space of quadratic differentials on (DS, Du) that 
are invariant with respect to the involution. We are most interested in a further subspace. 
For a Riemann surface (E,a;) with boundary, Q+(S,a;) is the subset of Q®(S,a;) that is 
non-negative everywhere on each boundary component, or equivalently where the boundary 
is horizontal (rather than vertical). For a surface S with no conformal structure, Q^(S) is 
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Data specified 



Length/Width 

Length 

Width 

Length and width 

Geometric 

Surfaces 

Gonformal 

structure 

MF+(E) 

MF-(E) 

Q+(E) 

Gombinatorial 

Graphs 

Elastic graph 

Length graph 

Width graph 

Strip graph 


Table 1. The different structures on surfaces and graphs 


the space of pairs of a conformal structure a; on E and a quadratic differential in Q+(E,a;), 
considered up to isotopy. Equivalently, a point in Q^(S) is a half-turn surface structure on S 
with horizontal boundary. 

From a quadratic differential q e Q+(E), we can get horizontal and vertical measured 
foliations 






At a point x e S, the vectors v e T^-S that are tangent to one of these measured foliations 
are those for which q{v,v) ^ 0 for q^ or q{v,v) ^ 0 for q"^. The transverse measure of an arc 
■y(t) is given by 


■= J|lmVg(7'W,YW)|c?^ ■= J|Re Vg(yW,yW) 


dt. 


A quadratic differential is more or less the combination of two of the other three types of 
data. 

• A conformal structure and a measured foliation T+ e AdJ-'^(E) uniquely determines 
a quadratic differential in g e with = q^. This is the Heights Theorem 
[HM79, KerSO, Gar84, MS84, Gar87]. 

• Let be a pair of measured foliations in AiF^(T,) and A4J-'^(S), respectively. 

Then generically there is a unique half-turn surface structure q e Q+(S) with the 
given foliations as horizontal and vertical foliations, respectively [GM91, Theorem 
3,1], 

• On the other hand, given a conformal structure and measured foliation T e A4T” (S), 
there is not always a quadratic differential in Q+ with F~ as its vertical measured 
foliation. You can always double the situation and consider the foliation D{F~) on 
i7(S), which has an associated quadratic differential restricting to a unique q e Q®(S), 
but q need not be in Q+(E); portions of dS might be vertical rather than horizontal. 
However, if there is such a quadratic differential, it is unique, as the doubling argument 
shows. 

The types of structure on surfaces are summarized on the top row of Table 1. 


3.2. Convexity. One key fact, used in Section 5.1, is that there is a natural convex structure 
on Recall that there are (several) natural coordinates for measured foliations. 

• On of a surface with boundary, pick a maximal collection of non-parallel 

disjoint simple arcs on E, and measure the transverse lengths of each. 
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• On of a closed surface or of a surface with boundary, take Dehn- 

Thurston coordinates with respect to some marked pair-of-pants decomposition of E. 
Normalize the Dehn-Thurston twist parameter so that twist 0 corresponds to measured 
foliations that are invariant under reversing the orientation of E. (See, e.g., [Thu].) 

Call any of these coordinate systems canonical coordinates. 

Definition 3.4. A function on Af J^^(E) is strongly convex if it is convex as a function on 
for each choice of canonical coordinates. 


This dehnition appears quite restrictive, since there are inhnitely many different canonical 
coordinates. However, such functions do exist. 


Theorem 2. For b e { + ,— }, let f be a function on weighted curves of type b so that 
• / does not increase under smoothing of essential crossings: 


(3.1) 




1 


where the crossing is essential and the strands have the same weight, and 
• f is convex under union: for C\ and C 2 two curves, 

(3.2) /[(C, u C2)/2] « (/[Cl] + f[C2])/2. 

Then f extends uniquely to a continuous, strongly convex function on AIJ-'^(E). 


An essential crossing of a curve is (somewhat loosely) one that cannot be removed by 
homotopy. Note that Ci u C 2 need not be simple, even if Ci and C 2 are. 


Example 3.5. Take any geodesic metric g on E, and let be the minimal length of 

a curve in the homotopy class [C] with respect to the metric g. Then Equation (3.1) is 
true, as smoothing crossings in the geodesic representative can only decrease length, and 
Equation (3.2) is true by dehnition (with equality). Thus all of these functions extend to 
continuous functions on AIJ-'^(E). Special cases of interest include when 

• the metric is a hyperbolic metric on E, or 

• the metric degenerates so that the lengths approach the transverse measure with 
respect to a measured foliation. 

As an example of this degeneration, the function F 1 —>■ i{\C^,F) for a hxed background 
curve C* is a strongly convex function on Af J-'^(E). 


Lemmas 5.4 and 5.5 below say that extremal length also satishes the same conditions. 

Proof sketch of Theorem 2. Fix a canonical set of coordinates on AIJ-'^(E). Given two rational 
measured foliations Fq, Fi e A1let Fi/2 = {Fq® Fi)/2 be the midpoint of the straight-line 
path between them, where © means adding the chosen canonical coordinates. (This depends 
on the coordinate system.) Since they are rational, Fq and Fi are represented by weighted 
simple curves, as is A 1 / 2 . Analysis of the coordinates shows that A 1/2 is obtained from 
(Fq u Fi)/ 2 by smoothing crossings (for some choice of resolution of the crossings). Thus, by 
Equations (3.1) and (3.2), 

f[Fy2] < f[{Fo U A)/2] < (/[Fo] + /[Fi])/2, 

which implies that /, when dehned, is a convex function. But / is dehned on the dyadic 
rational points in AfF^(E), and continuity of / follows from convexity [Roc70, Theorem 
10 . 1 ]. □ 
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Figure 3. 


Thickening a graph into a surface, T NT 



Figure 4. Geometrically thickening a graph. Left: An edge of length a is 
thickened to an a x e rectangle. Right: At a vertex, glue each half of the end 
of each rectangle to one of the neighbors. This is a more geometrically precise 
version of Figure 3. 


3.3. Graphs. In this paper, a graph F is a connected 1-complex, with possibly multiple 
edges and self-loops. A ribbon graph has in addition a cyclic ordering on the ends incident to 
each vertex; this gives a canonical thickening of F into an oriented surface with boundary 
A^F, as in Figure 3. 

A curve on a graph F is a map from a 1-manifold C into F, considered up to homotopy. 
We do not admit arcs here, and up to homotopy we can assume that curves are taut: they do 
not backtrack on themselves. 

We can put geometric structures on graphs corresponding loosely to the four geometric 
structures on surfaces above, as outlined in Table 1. 

First, corresponding to a conformal structure, we can consider an elastic graph T, with an 
(elastic) weight on the edges: for each edge e of F, give a positive measure which is absolutely 
continuous with respect to the Lebesgue measure. Up to equivalence, we effectively just give 
a total measure a(e) on each edge. 

There are at least two ways to interpret this measure a. On one hand, we can create a 
rubber-band network. For each edge e, take an idealized rubber band with spring constant 
1/Q;(e), so the Hooke’s Law energy when the edge is stretched to length ^(e) is 

(3.3) Energy(e) = t'(e)^/a(e). 

Attach these rubber bands at the vertices. Note that, as for real rubber bands, a longer 
section of rubber band (larger a{e)) is easier to stretch. Unlike for real rubber bands, the 
resting length is 0. 

Alternately, we can think of a{e) as dehning a family of rectangle surfaces: given elastic 
weights on a ribbon graph F and a constant e > 0, define a conformal surface NJ? by taking 
a rectangle of size a(e) x e for each edge e and gluing them according to the ribbon structure 
at the vertices, as indicated in Figure 4. T should be thought of as the limit of the Riemann 
surfaces N^T as e —> 0. 






FROM RUBBER BANDS TO RATIONAL MAPS 


17 


The next geometric structure is a width graph, which is a graph C with a width on the 
edges: for each edge e, there is a width w(e) satisfying a triangle inequality at each vertex v, 
as follows. Let Ci,..., e„ be the edges incident to v (with an edge appearing twice if it is a 
self-loop). Then we require, for each 1 < j < n, 

(3.4) w(ej) < 2 w(ei). 

l^i^n 

i^j 

W(r) is the space of possible widths on an abstract graph L. 

For a ribbon graph T, there is a natural surjective map W : (NT) W(r). For e 

an edge of F, let e* be the arc on A^F dual to e. Then W{F^) takes the value F^{e*) on the 
edge e. When F is trivalent (i.e., the vertices of F have degree < 3), fF is an isomorphism. 
In general there is ambiguity about how to glue at the vertices. 

Dually, a length graph K is a. graph together with an assignment i of a non-negative length 
to each edge of F, which we think of as defining a pseudo-metric on F. (It is a pseudo-metric 
because points can be distance zero apart.) T(F) is the space of lengths on F. 

If F is a ribbon graph with no valence 1 vertices, there is a natural injective map L* : 
/:(F) MJ^-(iVF), defined by 

eGEdge(r) 

for i a loop on F. This map is not surjective, but if F is trivalent, L* maps onto an open 
subset of AIJ-'^(iVF). 

Finally, a strip graph S' is a graph in which each edge e is equipped with 

• a length £(e), 

• a width w(e), and 

• an elastic weight (aspect ratio) a(e), 

so that i[e) = a{e)w{e) and the widths satisfy the triangle inequalities. Equation (3.4). A 
strip graph has a natural underlying elastic graph r(S'), width graph C'(S'), and length graph 
K{S). It also has a total area 

(3.5) Area(F) = ^ ^{e)w{e) = ^ ^ a{e)w{ef. 

eGEdge(S') eGEdge(5) ^ ' eGEdge(5) 

(Compare this to elastic energy. Equation (3.3) above, and extremal length on graphs. 
Equation (5.5) below.) 

Warning 3.6. Be careful to distinguish between lengths and elastic weights. They are 
determined by the same data (a measure or equivalently a length on each edge), but they 
are interpreted differently. If the elastic weights are interpreted as dehning a rubber band 
network, then the lengths can be interpreted as lengths of a system of pipes through which 
the rubber bands can be stretched. Alternatively, when the elastic weights are aspect ratios 
of rectangles (length/width), the lengths are just the length of the rectangles. 

As for surfaces, two of the three other types of data determine a strip graph, except that a 
choice of elastic weights and lengths may not correspond to a strip graph, as the triangle 
inequalities on the widths may be violated. 

Remark 3.7. In the thickening of an elastic graph F into a surface NJf in Figure 4, the 
precise details of how you glue the rectangles at the vertices are irrelevant in the limit, in 
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the following sense. If we pick two different ways of gluing at a vertex (e.g., gluing different 
proportions to the left and right) and get two families of surfaces N^T and N^T, then in the 
limit as e 0 the minimal quasi-conformal constant of maps between N^F and N^F goes 
to 1. 

In particular, a strip graph S has a thickening N^S, where each edge e is replaced by a 
rectangle of size i{e) x ew{e), and the rectangles are glued in a width-preserving way at the 
vertices. (This gluing is possible by the triangle inequalities.) If F{S) is the underlying elastic 
graph of S, then for e « I the surfaces N^S and Nf^F{S) are nearly conformally equivalent. 
See Proposition 5.10 for some related estimates. 

3.4. Extensions. The graphs T we are considering in this paper are usually topological 
spines of the corresponding surface NF (i.e., NF deformation retracts onto T). More generally, 
we may consider a graph T and surface E with a vri-surjective embedding T S, or more 
generally still a graph T, group G, and surjective map vri(r) ^ G (i.e., a generating graph 
for G). Much of the theory extends to this case. 

It is also sometimes convenient to generalize from surfaces to orbifolds. In the setting of 
groups, this means considering maps from 7ri(r) to the orbifold fundamental group 7r°’’'’(E). 
In the setting of graphs embedded in surfaces, we consider graphs with marked points that 
are required to be mapped to marked (orbifold) points. In particular, this is required for a 
proper statement of Theorem 1 if we drop the restriction that there be a branch point in 
each cycle of P. 


4. Harmonic maps and Dirichlet energy 
We now turn to harmonic maps, maps that minimize some form of Dirichlet energy. 

4.1. Harmonic maps from surfaces. Given a conformal Riemann surface (S,a;), a length 
graph and a Lipschitz map f : T, ^ K, the Dirichlet energy of / is 

(4.1) Dir(/):= f \Vf{z)\UA. 

Jt. 

Here, we have picked an (arbitrary) Riemannian metric g in the given conformal class, dA is 
the area measure with respect to g, and \Vf{z)\ is dehned to be the best Lipschitz constant 
of / at 2 ; with respect to g and i. (This agrees with the usual norm of the gradient when 
f{z) is in the interior of an edge of K and / is differentiable.) 

We also dehne Dir[/] to be the lowest energy in the homotopy class of /: 

Dir[/] := inf Dir(c/). 

9S[/] 

This optimum is achieved [EFOl, Theorem 11.1]. The optimizing functions are called harmonic 
maps and dehne canonical quadratic differentials in Q^(S). 

Warning 4.1. The connection between harmonic maps and quadratic differentials is not as 
tight as you might expect. Given a closed surface S and a weighted simple curve (7, we can 
construct a length graph K{C), with vertices the connected components of E\G and edges 
given by components of G, with lengths of edges of K{C) given by weights on G. There is a 
natural homotopy class of maps S ^ -^(G). The Dirichlet minimizer in this homotopy class 
sometimes, but not always, recovers the Jenkins-Strebel quadratic differential with vertical 
foliation given by G. 
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A similar construction using maps from the universal cover of S to an M-tree does recover 
the quadratic differential with given vertical foliation. Wolf used this to give an alternate 
proof of the Heights Theorem [Wol98]. 

Alternatively, we could generalize the target space, and allow A' to be a general non- 
positively curved polyhedral complex. In this context, we can take a complex K^{C) made 
of very thin tubes around the edges around the graph K{C) defined above. Then there is 
a closer connection between harmonic maps to K^{C) and the Jenkins-Strebel quadratic 
differential. We do not pursue this further here. 

4.2. Harmonic maps from graphs. Given an elastic graph (T, a) and a length graph 
{K,£), the Dirichlet energy of a Lipschitz map / : T ^ A' is 

(4.2) Dir(/) := f \f{x)\^dx= \ Fi\\f{y)dy 

Jr Jk 

(4.3) Dir[/] := irif Dir(^). 

9e[f] 

Here, \f'{x) \ is the derivative of / with respect to the natural coordinates given by a on T and 
£ on K. This derivative is not defined at vertices, but these points are negligible. Alternatively, 
\f'{x)\ is the best Lipschitz constant of / at x. Fill/(?/) is the filling function of f ai y e K, 
the sum of derivatives at preimages: 

(4-4) Fill/fe) 2 \f(x)\. 

f{x)=y 

The two expressions for Dir(/) are related by an easy change of variables. 

Minimizers of Dirichlet energy are harmonic functions in the following sense. 

Definition 4.2. A function / : T —> A' from an elastic graph T to a length graph K is 
harmonic if the following conditions are satisfied. 

(1) The map / is piecewise linear. 

(2) The map / does not backtrack (i.e., / is locally injective on each edge of T). 

(3) The derivative |/'(a;)| is constant on the edges of T when defined. As a result, for e 
an edge of T we may write \f{e)\ for the common value at any point on the edge. 

(4) If a vertex u of T maps to the interior of an edge of K, with edges ei,..., incident 
on the left and edges e^+i,..., e^+i incident on the right, then 

(4-5) 2l/'(ei)| = 2l/'(ei+fe)|. 

In particular, for each edge E of K, the filling function Fillj(r/) is constant on E, and 
we may write Fill/-(A). 

(5) If a vertex u of F maps to a vertex w of K, let Ai,..., A^ be the germs of edges of K 
incident to w. Then we have a vertex balancing condition: for i = 1, ..., k, 

(4.6) 2 V'{e)\ < T li 

e incident to v l^j^k e^-^Ej 

ei-^Ei j^i 

Here, the notation ” means the sum over all germs of edges of F that locally 

map to Aj. 

Theorem 4.3. A function f : T ^ K is a local minimum for the Dirichlet energy within its 
homotopy class iff it is harmonic. Every local minimum is also a global minimum. 
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Figure 5. A simple marked elastic graph and the associated tiling by rect¬ 
angles. All elastic weights on the graph are 1 except for one edge, which has 
weight 2. On the rectangle tiling, we have shown the widths of the rectangles, 
which are \f'{e)\ in the associated harmonic map to the interval. 


From the rubber bands point of view, the intuition behind Definition 4.2 is that 2\f'{x)\ is 
the tension in the rubber band, which is constant along the edge. The net force on a vertex 
mapping to the interior of an edge must be zero; this is condition (4). For a vertex v mapping 
to a vertex w, the net force pulling v in any one direction away from w cannot be too large; 
this is condition (5). Condition (4) can be thonght of as the special case of condition (5) 
when w has only two incident edges. 

Conditions (4) and (5) also imply that the derivatives \f'{e)\ form a valid width structure 
on r and that the hlling functions Fillj(F) form a width structure on K. From the rectangular 
surface point of view, a harmonic map gives a tiling of K with rectangles, with each edge e 
of F contribnting a rectangle of aspect ratio a{e). To illnstrate this, consider the case of 
marked planar graphs. Let F be a planar graph with two distingnished vertices s and t 
bordering the inhnite face, and consider the minimizer of the Dirichlet energy among maps 

(r,s,()^([o,i],{o),{i)) 

mapping F to the interval, taking s to 0 and f to 1. Brooks, Smith, Stone, and Tutte showed 
that this minimizer gives a rectangle packing of a rectangle [BSST40]. See Fignre 5 for an 
example. 

Remark 4.4. In the context of rectangle packings, it is traditional to nse the langnage of 
resistor networks rather than elastic or rubber-band networks. The Dirichlet energy in these 
two settings is the same for maps to an interval, but for general target graphs rubber bands are 
more flexible, because of the lack of orientations. Where rubber band optimizers are related to 
homotopy of maps from T to a target graph and holomorphic quadratic differentials, electrical 
cnrrent or resistor network optimizers are related to homology of F and to holomorphic 
differentials (not quadratic differentials). 
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4.3. Relation to Lipschitz energy. Another natural norm for maps cj): Ki ^ K 2 between 
metric graphs is the Lipschitz energy 

(4.7) Lip(0) := esssup|0'(x)| 

xgKi 

(4.8) Lip[0] := inf Lip(V>). 

he[0] 

Compared to Equation (4.2), we take the norm of |0Xa;)| rather than the norm. 

Let us note that 

• Lip(0) < 1 iff 0 is distance-decreasing, and 

• the Lipschitz energy is dynamical, in the sense that it behaves well under composition: 

Lip[0 o -0] ^ Lip[0] Lip['0]. 

The Lipschitz energy should be thought of as an invariant of a map between two length 
graphs, not as an invariant of a map from an elastic graph to a length graph. 

It follows immediately from the dehnitions that Lip is sub-multiplicative with respect to 
both length and Dirichlet energy: for C a curve on Ki and f : T ^ Ki a. map from an elastic 
graph, 

iK2[(l^*C] < Up[(j)]iKdC] 

Dir[0o/] ^ Lip[0]^Dir[/]. 

These are both special cases of Equation (1.8), and in both cases the inequalities are tight. 
For instance, we have [FMll, Proposition 3.11] 

.9) Lip[0] = sup 

C curve on Ki ^Ki 

4 . 4 . Computing Dirichlet energy. Given a homotopy class [/] : P ^ 77 of maps between 
an elastic graph P and a metric graph K, how can one hnd a representative g e [/] that 
minimizes Dirichlet energy? By Theorem 4.3, this is the same as Ending a harmonic map in 
the homotopy class. 

As with harmonic maps in other settings, this is easy to compute, with at least two 
reasonable approaches: 

(1) Repeated averaging, iteratively moving the image of each vertex of P to the weighted 
average of its neighbors in the universal cover K of K. The average of a non-empty 
set S of points on a metric tree is the (unique) point x that minimizes the sum of 
squares of distances from x to points in S. Note that this average can generically be 
on a vertex. 

(2) Linear or convex quadratic programming, based on the observation that once the 
combinatorics of the map are hxed, specifying which vertices of P go to which vertices 
or edges of K, the Dirichlet energy is a convex, quadratic function of the positions 
along the edges, and thus the minimum can be found by solving linear equations. 
Here, one starts by guessing some combinatorics, and updating the combinatorics if it 
turns out not to be optimal. 

The second approach should usually be faster, but also requires some more care in changing 
the combinatorics. The harmonic representative of [/] is not unique in general; however, the 
set of harmonic representatives forms a convex set in a suitable sense. 
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5. Extremal length 


5.1. Extremal length on surfaces. Given a Riemann surface (E,a;) and a curve C on S, 
recall that the extremal length of G on S is (omitting analytic details) 


(5.1) 


EL(^[G] := sup 


p:S- 


^0 Areapg(E)’ 


where we make the following definitions. 

• The metric g is an arbitrary metric in the conformal class u. 

• The metric pg is the metric g scaled by the conformal factor p. (This may be a 
pseudo-metric.) 

• The number ^pp[G] is the minimal length of any element of the homotopy class [G] 
in the metric pg. 

• The number Areapg(S) is the total area of S with respect to pg. 

Observe that scaling p by a global constant does not change the supremand. 

We can extend Equation (5.1) to allow for weighted curves: for G = ^ weighted 

curve, define £[G] := use Equation (5.1) as before. (Note that this agrees 

with the earlier definition for integral weights.) 

The following results are standard. 


Theorem 5.1 (Jenkins-Strebel). If C is a weighted simple curve, then the supremum in the 
definition o/EL[G] is achieved, and the optimal metric is the metric on a half-turn surface 
associated to a guadratic differential with horizontal foliation equal to C. 

Lemma 5.2. Let vr : E ^ E be a covering map of degree d. For C a weighted curve on E, 
define 7r~^C to be the inverse image of C, with the same weights. Then EL[7r“^G] = (iEL[G]. 

Lemma 5.3. If a e M>o is a global weighting factor, then 

EL[aC] = a^EL[C]. 

We now turn to properties related to convexity, as in Section 3.2. The following two 
properties follow from elementary arguments. 


Lemma 5.4. EL does not increase under smoothing of essential crossings: 



Here, the two sides show a local picture of unweighted curves, or weighted curves with equal 
weights on the two local strands. 


Lemma 5.5. If Ci and C 2 are two weighted multi-curves, then 

EL[Gi u G 2 ] ^ 2(EL[Gi] + EL[G2]), 

where we keep the weights on each component of Ci and G 2 . More generally, for 0 < t ^ 1, 

EL[tGi u (1 - t)C2] < tEL[Gi] + (1 - t) EL[G2]. 

Remark 5.6. In Lemma 5.5, we are taking the union of curves, not the union of path families, 
as sometimes appears in the theory of extremal length. 

Corollary 5.7. EL extends uniquely to a continuous, strongly convex function on A1J^*(E). 
Proof. Follows from Theorem 2. □ 
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Remark 5.8. Presumably the techniques in the proof of Corollary 5.7 can be extended to 
prove the rest of the Heights Theorem, in particular that the associated quadratic differential 
varies continuously as a function of the measured foliation. 


5.2. Extremal length on graphs. By analogy with Equation (5.1), for C a weighted curve 
on an elastic graph (T, a), dehne 


(5.2) 


EL„[C] 


sup 

p:Edge(r)^]R5.0 


Areap„(r)’ 


where 


• The length metric pa on T gives edge e the length p(e)a(e). 

• The number lpa\C] is the length of C with respect to pa, i.e., 

(5.3) ^po\C] := Yj ^c(e)p(e)a(e), 

eeEdge(r) 

where nc{e) is the weighted number of times that C runs over e. 

• AreapQ,(r) is the “area” of T with respect to pa, defined to be 

(5.4) Areap«(r) := Y p(e)^a(e). 

eGEdge(r) 

(The intuition is that each edge is turned into a rectangle of width proportional to 
p(e) and aspect ratio a(e), and thus area p(e)a(e).) 

In fact the supremum in Equation (5.2) is easy to do. The optimum has p(e) proportional 
to nc(e) and so 


(5.5) 


EL„[C] = 2 nc{efa{e). 

eGEdge(r) 


This formula extends immediately to a function on VV(r), and satishes Lemmas 5.2, 5.3, 
and 5.5. Eor Lemma 5.4, we need to pick a ribbon structure on T in order to make sense of 
“essential” crossings. With any such choice. Lemma 5.4 is true. 


Corollary 5.9. Let (T, a) be an elastic spine for a surface S. Then the function C ELq,[C'] 
on C^(S) extends uniquely to a strongly convex function on AfJ^^(S). 


5.3. Relating graphs and surfaces. We can now give a concrete relation between an 
elastic graph (T, a) and the associated family of conformal surfaces NJ?. Write EL[C; T] for 
extremal length with respect to the elastic graph T, and EL[C'; E] for extremal length with 
respect to the conformal surface E. 


Proposition 5.10. Let (T, a) be an elastic ribbon graph with trivalent vertices, and let 
m = minea(e) he the smallest weight of any edge in T. Then, for t < m/2 and C any 
measured foliation onT, we have 

EL[C;r] <fEL[C';Wr] ^ EL[C'; T] • (1 + 8f/m). 

The proof involves finding, on the one hand, embeddings of sufficiently thick annuli into 
AfT, and, on the other hand, suitable test functions p on NtV in Equation (5.1). 


Remark 5.11. The restriction to trivalent graphs in Proposition 5.10 can be removed. Note 
that the estimate depends only on the local geometry of T, and thus is unchanged under 


covers. 
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5.4. Duality with Dirichlet energy. As mentioned earlier, extremal length is in some 
sense dual to Dirichlet energy. More precisely, we have the following. 

Proposition 5.12 (Sub-multiplicative). Let C be a curve in an elastic graph F, and let 
f \ T ^ K he a harmonic map to a length graph. Then 

i[f,Cf ^I)iT[f] EL[C]. 

Proof sketch. This is basically the dehnition of EL in Equation (5.2), noting the similarities 
between the area in Equation (5.4) and Dirichlet energy. Equation (4.2). □ 

Proposition 5.13 (Duality 1). Let f \ T ^ K he a harmonic map from an elastic graph to 
a length graph. Then there is a sequence of weighted curves Ci in F so that for all e > 0 there 
is an i so that 

Dir[/] EL[Q]. 

Proof sketch. Take weighted curves Ci so that ncXe) approximates |/'(e)|. □ 

Proposition 5.14 (Duality 2). Let C be a curve in an elastic graph F. Then there is a 
length graph K and a harmonic map f : T K so that 

^{f,Cf = Biv{f] EL[q. 

Proof sketch. Take A to be F with edge lengths a{e)nc{e). □ 

The situation is less satisfactory for surfaces. Sub-multiplicativity in the sense of Propo¬ 
sition 5.12 is true (at least when the curve C is embedded), as is Proposition 5.13. Issues 
related to Warning 4.1 make an analogue of Proposition 5.14 more delicate, although it is 
true if we allow non-positively curved polyhedral complexes as the target space. 

Remark 5.15. The dehnition of extremal length on graphs in Equation (5.2) is somewhat 
backwards, in that it is a supremum of a ratio of energies over all metrics (or equivalently 
over all maps to length graphs). By analogy with Dirichlet energy (Equation (4.2)), it would 
be better to dehne the energy of a homotopy class as an inhmum of some energy functional. 
Indeed, we could take Equation (5.5) as the primary dehnition. 

This remark applies to extremal length on surfaces (Equation (5.1)) as well: it might be 
better to take a different dehnition as primary. Namely, recall that a conformal annulus A 
has an extremal length EL(A), which we may dehne as the inverse of the modulus. Then 
extremal length of a weighted multi-curve C = Yji O'iCi can be alternately dehned as 

(5.6) EL[C'] = inf V EL(Ai), 

I 

where the inhmum runs over all disjoint embeddings of conformal annuli Ai with core curves 
homotopic to C*. 

5.5. Extensions. There are several ways in which we can extend these notions of extremal 
length. Eirst, we can consider graphs that are embedded in a surface, not necessarily as a 
spine. 

Definition 5.16. A graph F embedded in a surface S is filling if each component of S\F is 
a disk or an annulus on the boundary of S. 
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Definition 5.17. Let {K,i) a length graph with a filling embedding 0 in a snrface E. For 
C a curve in S, define 

£^[C']:= inf J[C']. 

C curve on K 

Now for (r,a) a filling elastic graph in S, define ELr,o[C'] by 

^r,pa[Cf 

A-rbV- 

p:Edge(r)—/ 

This is just like Equation (5.2), except that we consider homotopy classes in S rather than 
in r. A similar notion of extremal length was considered by Duffin [Duf62], in the context of 
electrical networks and graphs with two marked points (as in Section 4.2). 

The optimization in Definition 5.17 is no longer as easy, and the analogue of Equation (5.5) 
is more awkward to state. The result of the optimization gives a rectangular tiling of the 
surface with aspect ratios given by a, analogous to Figure 5. 

More generally, we can consider graphs generating a group. 


ELr,^^] : 


Definition 5.18. Eor (F, a) an elastic graph, 0 : 7ri(F) 
onto a group G, and {g\ a conjugacy class in G, define 

inf J-pa[C] 

urve on F 

<^[C] = [9] 


%a\3\ ■ „ 

C curve on i 


G a surjective homomorphism 


EL'Jl?] 

EL„[ 5 (] 


sup 

p:Edge(r)—>IRj.o -^i'eapQ,(F) 
hm ELdl7"]/n'. 

n^oo 


There are also other models for defining a combinatorial length of curves on a graph. 
Notably, Schramm [Sch93] and Canon, Floyd, and Parry [CFP94] define a model where the 
length of a curve in a graph is determined by the vertices that it passes through, rather than 
the edges that it crosses over (as in this paper). More generally, one can consider shinglings 
of a graph or surface, decompositions of the space into a finite number of overlapping open 
sets. The combinatorial length of a curve is then determined by which shingles it passes 
through. The edge model that is the main focus of this paper comes from taking shingles 
that are neighborhoods of the edges (overlapping at the vertices), while the vertex model of 
Schramm-Cannon-Floyd-Perry comes from taking shingles that are neighborhoods of the 
vertices (overlapping at the centers of edges). 

For any of these notions of length, one can define a notion of extremal length using 
Equation (5.2). 

More generally, one may instead attempt to characterize which extremal length functions 
can appear. There are several natural sources of an “extremal length” function on AdJ^^(S): 

(1) A conformal structure a; on S gives the usual notion of extremal length. 

(2) If E' is another surface and 0 : E' E is a filling embedding (an embedding for 
which the complementary regions are disks or annuli), then a conformal structure u' 
on E' gives a notion of extremal length on E, defined analogously to Definition 5.17. 

(3) An elastic graph spine (F, a) for E gives a notion of extremal length by Equations (5.2) 
or (5.5). 

(4) A filling elastic graph in E gives a notion of extremal length by Definition 5.17. 

(5) Finally, a shingling of E gives yet another notion of extremal length. 
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All of these notions of extremal length give a function on C+(S) that 

• is positive, 

• is homogeneous quadratic, in the sense of Lemma 5.3, 

• does not increase under smoothing, in the sense of Lemma 5.4, 

• is sub-additive under union, in the sense of Lemma 5.5, and therefore 

• extends to a strongly convex function on by Theorem 2. 

As a result of convexity, we can think of EL as a kind of “norm” on 

Problem 5.19. Which functions EL : ^ M can arise from the constructions 

above? 

The properties above are some restrictions, but are probably not a complete list. 

Some of these notions of extremal length subsume the others: Extremal lengths from notion 
(5) include extremal lengths from notion (4) by a direct construction. By Proposition 5.10, 
notion (2) is dense in notion (4) (up to scale). Notion (4) naturally includes notion (3), and 
by taking the graph to be the edges of a triangulation we can see that notion (4) is dense in 
notions (1) and (2). 

Problem 5.20. The dehnition of extremal length in surfaces. Equation (5.1), extends to 
(width) graphs embedded in E, rather than just curves. What do the resulting optimal 
metrics look like? 


6. Stretch factors and embedding energy 


Now we study how extremal length and Dirichlet energy change under maps between 
surfaces and graphs. This material is developed more fully in a sequel paper with Pilgrim 
and Kahn [KPT15]. 

6.1. Stretch factors for surfaces. Let 0 : Si ^ S 2 be a topological embedding of surfaces. 
Then composing with (j) and deleting null-homotopic components induces a natural map 
0* : C+(Ei) —> C+(E 2 ). (This pushforward does not work on C~(Ei).) 

Warning 6.1. The map (j)^ does not generally extend to a continuous map from AlJ^^(Si) to 
M.T+(E2). 

Now suppose El and E 2 have conformal structures Ui and a; 2 , respectively. (The map 0 
need not respect the conformal structures.) 


Definition 6.2. In the above setting, the stretch factor of (f is 


( 6 . 1 ) 


SF[0] := sup 

CsC+(E) 


This depends only on the homotopy class of <p. 


EL^JC-] 


It follows from the dehnition that SF behaves well under composition. 


Proposition 6.3. // / : Ei E 2 and : E 2 E 3 are two topological embeddings of 
conformal surfaces, then 

SF[fog]^SF[f]-SF[g]. 

Definition 6.4. A conformal embedding 0 : Ei ^ E 2 is strict if, in each component of E 2 , 
there is a non-empty open subset in the complement of the image. 
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Theorem 3 (essentially Ioffe). If (Ei,a;i) and (T, 2 ,(^ 2 ) are two conformal surfaces and 
0 : El E 2 is a homeomorphism so that there is no strict conformal embedding in the 
homotopy class [0], let K he the lowest constant so that there is a K-quasi-conformal map 
in the homotopy class [0]. Then 

SF[0] = K. 

Proof. This is very close to the main theorem of [Iof75]. In that paper, Ioffe proves that if 
there is no conformal embedding in [ 0 ], there are canonical quadratic differentials qi e Q+(Sj) 
and a iF-quasi-conformal representative for 0 that uniformly stretches qi to ^ 2 - Suitably 
approximating the horizontal foliations of Qi and Q 2 by rational measured foliations (being 
careful about Warning 6.1) gives a sequence of simple curves (not necessarily connected) so 
the ratio of extremal lengths approaches K. 

The case when there is a conformal embedding, but not a strict conformal embedding, can 
be treated by, for instance, adding annuli to the boundary components of Ei so there is no 
conformal embedding. □ 

Corollary 6.5. //Si and S 2 are closed surfaces and 0 : Si —> S 2 is a homeomorphism, then 
the Teichmuller distance between Si and S 2 is /logSF[0], in the sense that if Tq is a fixed 
base surface and 0 : Eq ^ Ei is a marking, then the distance between the marked surfaces 
(Si,0) and (S 2 , 0 o 0 ) is /logSF[0]. 

Proof. Immediate from Theorem 3 and the definition of Teichmuller distance. □ 

In this context. Proposition 6.3 is the triangle inequality for Teichmuller distance. Corol¬ 
lary 6.5 was proved by Kerckhoff [KerSO, Theorem 4]. In this case connected simple curves 
suffice. 

Remark 6 . 6 . Despite Warning 6.1, there is a non-continuous extension of 0^)= to measured 
foliations. Dehne a pull-back function 0* ; C“(E 2 ) —> l^C^(Si), where 0^C“(Si) is the set 
of subsets of C“(Si), by taking all possible intersections of C with the image of Si: 

0*(C) := { [c u 0(Si)] I c a representative of C }, 

where we delete inadmissible components as in the dehnition of 0^)=. Then, for a measured 
foliation Fi e AlJ^^(Ei), we may dehne 0 *(F 2 ) to be the (unique) measured foliation 
F 2 e M.F^{fp 2 ) so that, for all C 2 e C“(S 2 ), 

^([C 2 ],F 2 )= inf z([Ci],Fi). 

C\^(p^{C2) 

With this dehnition, the curves in Equation (5.1) can be replaced by measured foliations, 
and the supremum is achieved. 

Definition 6.7. The embedding 0 is annular if it extends to an embedding of an annular 
extension Si in S 2 , where Si is obtained by attaching an annulus to each boundary component 
of Si. 

Theorem 4 (Joint with Kahn and Pilgrim [KPT15, Theorems 1 and 2]). //(Si,a;i) and 
(E 2 ,a; 2 ) are two conformal surfaces and 0 : Ei E 2 is a topological embedding, then 0 is 
homotopic to a conformal embedding iffSF[(j)] < 1. 

Furthermore, the following conditions are equivalent: 

(1) SF[0] < 1, 

(2) 0 is homotopic to an annular conformal embedding, and 
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(3) (p is homotopic to a strict conformal embedding. 


Proof sketch. The case when cp is not homotopic to an embedding is implied by Theorem 3. 
The converse of the hrst claim is essentially Schwarz’s Lemma. 

To show that (1) implies (2), pick a quadratic differential q e Q+(Ei) that is strictly 
positive on each boundary component. Define to be Ei plus an annulus of width t on 
each boundary component, using the coordinates from q. Elementary estimates show that 
SF[S5^ ^ El] (in the natural homotopy class) approaches 1 as t approaches 0. Then for t 
small, by Proposition 6.3 

SF[S* ^ Ea] < SF[S* ^ Ei] • SF[Ei ^ E 2 ] < 1, 


so by Theorem 3 there is a conformal embedding of E^ in E 2 . 

Clearly (2) implies (3). Finally, if 0 is a strict conformal embedding of Ei in E 2 , then by 
considering test metrics we can show that 


( 6 . 2 ) 


SFW < 


sup 

geS+lSa) 

g#0 


Areaq(lm(0)) 

Areag(E2) 


Here lm(0) c E 2 is the image of cp, which by hypothesis misses an open subset of E 2 , 
and Areaq is the area with respect to the quadratic differential q. Thus for each q, 
Areag(lm(0))/Areag(E 2 ) < 1. Since the supremand doesn’t change as we scale q, we are 
maximizing over the compact set PQ+(E 2 ,a; 2 ) and the supremum is strictly less than 1. □ 


6.2. Behaviour of stretch factor under covers. In Section 7.2, we will also need an 
understanding of the behavior of SF under covers. Let 0 : Ei E 2 be a topological 
embedding of conformal surfaces, let E 2 be a finite cover of E 2 , and let 0 : Ei E 2 be the 
pull-back of 0; that is, Ei is defined by the pull-back of E 2 and Ei via the diagram 


El^ ^ E2 

(6.3) ; 

El E2. 


When E 2 is connected, Ei is the minimal cover of Ei with a map to E 2 making the above 
diagram commute. 

Question 6.8. How does SF[0] compare to SF[0]? 

It appears that SF[0] ^ SF[0] in general, at least if we allow E 2 to be an orbifold. However, 
we can make many partial statements. 

Lemma 6.9. For 0 and 0 as above, SF[0] ^ SF[0]. 

Proof. Follows from the definition of SF and the good behavior of extremal length under 
covers. Lemma 5.2. □ 

Lemma 6.10. For 0 and 0 as above, z/SF[0] ^ 1, then SF[0] = SF[0]. 

Proof. Follows from Theorem 3. □ 

Lemma 6.11. For 0 and 0 as above, SF[0] < 1 iff SF[(p] < 1 and SF[0] = 1 Zj^SF[0] = 1. 
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Proof. If SF[0] < 1, by Theorem 4 the map 0 is homotopic to a strict conformal embedding, 
and so 0 is too and thus SF[0] < 1. The rest follows from Lemma 6.10. □ 

Definition 6.12. For 0 : Ei E 2 a topological embedding of surfaces, dehne the lifted 
stretch factor by 

SF[0] := ^ lim SF[0], 

(p covers cp 

where the limit runs over increasingly large hnite covers of 0 (determined by a covering of 
E 2 ). These covers form a directed system, and SF only increases in a cover (Lemma 6.9) 
while remaining bounded (Lemmas 6.10 and 6.11), so the limit exists. 

By dehnition, if 0 is a covering of 0, then SF[0] = SF[0]. 

We will ultimately need to show that certain maps have lifted stretch factor less than one. 
We give two arguments, one more elementary and one proving a stronger result. (See the two 
proofs of Proposition 7.8.) 

Definition 6.13. A homotopy class of topological embeddings [0] : Ei ^ E 2 is conformally 
loose if, for all y e E 2 , there is a conformal embedding 0 e [0] so that y f 0(Ei). 

If E 2 is compact and [0] is conformally loose, then we can hnd hnitely many conformal 
embeddings 0* e [0], i = 1,..., n so that 

n 

(6.4) P 0(^=0. 

Proposition 6.14. //[0] is conformally loose wih n maps (pi in Equation (6.4), then SF[0] < 

1 — 1/n. 

Proof sketch. Fix a weighted multi-curve C on Ei, and let q e Q+(E 2 ,a; 2 ) be the quadratic 
differential corresponding to 0(0). For at least one i, we will have 

Areag(00Ei)) ^ ^ ^ / 

Areag(D) ^ 

Choosing test metrics as in Equation (6.2) shows that EL^j[0] ^ (1 — 1/n) EL^^^ [0*0], as 
desired. □ 

Corollary 6.15. //[0] is conformally loose, then SF[0] < 1. 

Corollary 6.15 is enough to prove Theorem 1. However, more is true. 

Theorem 5 (Joint with Kahn and Pilgrim [KPT15, Theorem 3]). For every strict conformal 
embedding 0 : Ei E 2 , 

SF[0] < 1. 

By Equation (6.2), Theorem 5 is a consequence of the following proposition. 


Proposition 6.16. Let E 2 be a Riemann surface, and let Ei c E 2 6e a subsurface with 
compact closure. Then there is a constant K < 1 so that, for all 0 q e Q(E 2 ), 

(6.5) < A, 

Areag(zj2) 

Furthermore, K can be chosen so that Equation (6.5) holds for all coverings of the pair 
(Si,F2). 
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Figure 6. The schematic setup of Lemma 6.17. 

Proposition 6.16, in turn, depends on the following local lemma. 

Lemma 6.17. Let Q V) be an open subset of the disk with an open set A in the complement 
ofQ, and let B B) be a neighborhood of Q r\ dD. Then, for every e > 0, there is a 6 > 0 so 
that, if q ^ Q(D) is such that 

( 6 . 6 ) 

then 

(6.7) 

Essentially, Lemma 6.17 says that if the area of q is concentrating in a subset of the disk, 
then it is concentrating near the boundary. See Figure 6. 

Proof sketch for Lemma 6.1 7. If there are no such bounds, there is an e so that we can End 
a sequence of quadratic differentials ^ Q(D) so that 


(6.8) 

Areag„(D) = 1 

(6.9) 

Areag„(5) < e 

(6.10) 

Areag„(12) > 1 


Consider |g„| as a measure on D. Since the space of measures of unit area on the closed 
disk is compact in the weak topology, after passing to a subsequence we may assume that 
\qn\ converges in the weak topology to some limiting measure /i on D. Absolute values of 
holomorphic functions are closed in the weak topology, so the restriction of /i to the open disk 
can be written |goo| for some holomorphic quadratic differential goo. But Areag„(A) < 1/n, so 
Areag^(A) = 0, so goo is identically 0; hence fi is snpported on dD. Eqnation (6.10) implies 
that the snpport of p is also contained in 12, and hence in 12 n dD. Bnt this contradicts 
Equation (6.9). □ 

Proof sketch for Proposition 6.16. Divide E 2 by smooth arcs aj so that the complementary 
regions are all disks; let these disks be t/j. Consider a qnadratic differential on E 2 with a 
very large proportion of its area in Ei. Then most t/j (as weighted by Areaq(17i)) must have 
most of their g-area in 27j n Ei. Lemma 6.17 then says that most of the area of most of the 
Ui mnst be in a small neighborhood of the seams a^. Arrange the constants so that more 
than half the total area must be in these neighborhoods. 


Areaq(12) 

T-> 1 - d. 

Areag(D) 

Aied.q{B) 

-——^ > 1 — e, 

Areaq(D) 




FROM RUBBER BANDS TO RATIONAL MAPS 


31 


Now pick an alternate set of seams I3k with the same properties, but so that j3k n = 0 
for all j and k. By the same argument, more than half the total area must be concentrated 
in a neighborhood of the Pk, a contradiction. 

All of the estimates in this proof only depend on the local geometry of the f/j, and thus 
remain unchanged under taking covers. □ 


6.3. Stretch factors and embedding energy for graphs. We now turn to the (easier) 
parallel theory of stretch factors for maps between graphs. Let 0:ri^r2bea continuous 
map between two graphs, and suppose that Li and r 2 have elastic structures oi and 02 , 
respectively. 


Definition 6.18. In this setting, the EL stretch factor of 0 is 

EL„,[0,C'] 


( 6 . 11 ) 


SFelW := sup . 

C curve on Fp J 


As for surfaces, this only depends on the homotopy class of 0. 
Dually, the Dirichlet stretch factor of 0 is 

Dir[0 o /] 


( 6 . 12 ) 


SFDir[0] := sup 

K length graph J 


where the supremum runs over all length graphs {K, i) and all homotopy classes of maps 
from r 2 to iL. 


Remark 6.19. As mentioned in the introduction, the Dirichlet stretch factor has a natural 
interpretation in terms of rubber-band networks. If SFDir[0] < 1, then the rubber-band 
network Ei is “looser” than the rubber-band network r 2 : minimal Dirichlet energy in a 
homotopy class with target an arbitrary length graph decreases under composition with 0. 
Theorem 6 below implies that this is true for arbitrary target geodesic spaces, not just graphs. 

As in the case of surfaces, these stretch factors behave well under composition (Proposi¬ 
tion 6.3). Unlike in the case of surfaces, for graphs we have a direct characterization of the 
stretch factor. 


Definition 6.20. For 0 ; Fi —> r 2 a Lipschitz map between elastic graphs, the embedding 
energy of 0 is 

(6.13) Emb(0) := esssupFil4(?/) 

yeT2 

(6.14) Emb[0] := inf Emb(0), 

where Fill^ is the hlling function dehned in Equation (4.4). We take the norm of Fill,^ 
rather than the norm used for Dirichlet energy. This is equivalent to Equation (1.1). 

By comparing Equations (1.1) and (4.7), we see that 

(6.15) Emb(0) ^ Lip(0). 

Theorem 6. For 0 : Ti —+ r 2 a continuous map between elastic graphs, 

SFel[0] = SFDip[0] = Emb[0]. 
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Theorem 6 should be thought of as analogous to Theorem 3, although it applies in all 
cases, not just when there fails to be a conformal embedding. It is also analogous to the 
relation between Lipschitz energy as maximnm derivative (Equation (4.7)) and as ratio of 
curve lengths (Eqnation (4.9)). 

Proof sketch. For C a curve on Fi and 0 : Fi ^ F 2 a Lipschitz map, by pushing forward the 
scaling function p in Eqnation (5.2) we can see that 

(6.16) EL[0,C'] < EL[C'] • Emb(0). 

This immediately implies that SFel[(;/)] < Emb[(;/)]. 

Similarly, the dehnition of the Dirichlet energy as the norm of the hlling fnnction 
and the embedding energy as the norm of the hlling fu n ction makes it clear that, for 
f ■.V 2 ^ K any map to a length graph, 

(6.17) Dir(/ o (j)) ^ Dir(/) Emb(0), 
which implies that SFDii.[0] < Emb[(;/)]. 

To prove the opposite inequalities, we hnd and analyze a representative for the homotopy 
class [0] for which Equations (6.16) and (6.17) are tight. We defer the description of these 
nice representatives (analogous to harmonic maps, Dehnition 4.2) nntil Section 6.5. □ 

Motivated by Theorem 6, we make the following dehnition. 

Definition 6.21. For (Fi,ai) and (F 2 ,a 2 ) two elastic strip graphs, a map 0 : Fi ^ F 2 is 
loosening if Emb(0) < 1. The map is strictly loosening if Emb(0) < 1. We likewise say that 
a homotopy class [0] is (strictly) loosening if there is a (strictly) loosening map in [0]. 

In contrast with the case for surfaces (Section 6.2), it is easy to see how the stretch factor 
for graphs behaves nnder covers. 

Proposition 6.22. Let 0 ; Fi ^ F 2 he a map of elastic graphs, F 2 be a cover 0 /F 2 , and 
0 : Fi —> F 2 he the pull-back map, defined as in Eguation (6.3). Then 

SF[?] = SFW. 

Proof. The dehnition of SF[0] by the maximal expansion of extremal length or Dirichlet 
energy shows that SF[0] ^ SF[0] (as in Lemma 6.9). The dehnition of Emb[0] by the 
inhmum over all representatives shows that Emb[0] < Emb[0]. Theorem 6 shows that we 
have eqnalities. □ 

6.4. Relating graphs and surfaces. As motivation for the somewhat strange dehnition 
of embedding energy, consider two elastic graphs (Fi,ai) and (F 2 ,a 2 ) and a conformal 
embedding of the thickenings WEi A"tF 2 . Snppose that this conformal embedding is 
“close” in some sense to a graph map 0 : Fi —> F 2 . If we look away from the thickenings 
of the vertices (of both graphs), we see locally a map from some edges of Fi into a single 
edge of F 2 . The total height of images of the thickened edges of Fi mnst be less than or 
eqnal to the total available height in the thickened edge of F 2 . Near a point x e Fi, the 
image is stretched horizontally by a factor of |0'(x)|; since we are considering a conformal 
embedding, the image mnst also be stretched vertically by the same factor. Thus the image 
of this portion of the edge takes up a height of f|0'(x)|. This argument suggests that if there 
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Ti 


r. 


Thicken 


V/////////////A 


N.Fi 


1—> 



A^.r2 


Figure 7. Motivation for the definition of embedding energy. Two edges of 
Ti mapping to one edge of r 2 get thickened np to two rectangles mapping to 
one rectangle. For the result to be a conformal embedding, we need the total 
height of the image to be less than the total height available in the range. 

is a conformal embedding close to 0 , we must have, for each y ^T 2 not near a vertex or the 
image of a vertex, 

^ t or, equivalently, Fill 0 (i/) ^ 1 . 

See Figure 7. 

We will not attempt to make the above heuristic argument precise. Instead, we get a 
precise statement another way. 

Proposition 6.23. Let Fi and r 2 be trivalent elastic ribbon graphs, and let m = mineQ!(e) 
be the smallest weight of any edge in either Fi or r 2 . Let 0 :ri^r2 be a map that extends 
to a topological embedding Ntf): WFi —» NtT 2 . Then, for t < m/2, we have 

SF[ 0]/(1 + 8t/m) < SF[W0] ^ SF[ 0 ] • (1 + 8t/m). 

Proof. Immediate from Proposition 5.10. □ 

Theorem 7. Let Fi and T 2 be elastic ribbon graphs, and /ef 0 : Fi ^ F 2 be a map of graphs 
that extends to a topological embedding Ncf : iVFi —> iVF 2 . //Emb[0] < 1, then for all t 
sufficiently small, WFi conformally embeds in WF 2 in [-/V0]. Conversely, if NtVi conformally 
embeds in NtV 2 in [iV0] for all sufficiently small t, then Emb[0] < 1. 

Proof. Follows from Proposition 6.23 and Theorem 4. If the graphs Fi and F 2 are not trivalent, 
we may approximate them by trivalent graphs at the cost of a small change in the stretch 
factor. □ 

6.5. Minimizing embedding energy: A-filling maps. We say a little more about the 
proof of Theorem 6 , and in particular characterize the optimal maps. Recall from Section 3.3 
that a strip graph has both lengths and widths. 

Definition 6.24. For A > 0, a map (f : Si ^ S 2 between strip graphs is X-filling if 

(1) 0 is length-preserving: it does not backtrack, and 10 ^ 2 :) | = 1 for almost all x, where 
we take derivative with respect to the length metrics, and 
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(2) 0 scales widths by A: for almost all y e S 2 , 

Yj w{x) = A • w{y). 

Note that for maps between strip graphs, we have to be careful whether we differentiate 
with respect to the length metric or with respect to the coordinates from the elastic weights, 
as in Warning 3.6. 

Lemma 6.25. //0 : S'! —> S '2 is a X-filling map between strip graphs, then the underlying 
map r(0) ; r(S'i) —> r(S' 2 ) between elastic graphs has embedding energy A. This is minimal 
in the homotopy class. 

Proof sketch. The conditions on (j) imply that Fillr(<;i) = A everywhere on r(S' 2 ). (Recall that 
“length” on the elastic graph is length/width in terms of the strip graph.) 

Now consider the length graph K{S 2 ). The map r(S' 2 ) —> K{S 2 ) has Dirichlet energy 
Area(S' 2 ), almost by dehnition. The composite map r(S'i) —> r(S' 2 ) ^ K{S 2 ) has Dirichlet 
energy A ■ Area(S' 2 ). It follows that SFDir[r(0)] ^ A. The easy direction of Theorem 6 then 
implies that Emb[F(0)] = A. □ 

We say that a map (j) between elastic graphs is X-filling if there are compatible strip 
structures that make (j) A-hlling. It is not true that every homotopy class of maps between 
elastic graphs has a A-hlling representative. However, we can make it A-hlling on a subgraph. 

Definition 6.26. A map : Si S 2 between strip graphs is partially X-filling if there are 
non-empty subgraphs Ti of Si and T 2 of S 2 so that 

(1) 0(Ti)=T2and fi-\T2) = Tp, 

(2) the restriction of 0 to a map Ti T 2 is A-hlling; 

(3) 0 is everywhere length-preserving; and 

(4) outside of Ti and T 2 , the map cj) scales widths by less than A. 

Proposition 6.27. Every homotopy class of maps between elastic graphs has a partially 
X-filling representative. 

Theorem 6 follows quickly from Proposition 6.27. 

6.6. Computing embedding energy. Given two elastic graphs Fi and F 2 and a homotopy 
class [0] : Fi ^ F 2 of maps between them, how can we concretely hnd the partially A-hlling 
representative fi e [0] guaranteed by Proposition 6.27 and thus compute Emb[0]? The 
following iteration appears to converge. 

(0) Pick a set of widths vq e W(F 2 ), i.e., a width for each edge of F 2 satisfying the triangle 
inequality. 

(1) Take the metric graph Kq to be P 2 with edge e assigned length a{e)vo{e). Note that 
the evident map /o : F 2 Kq is harmonic. 

(2) Eind a harmonic representative go of the composite map Fi F 2 ^ Kq. Our hrst 
approximation to fi is fio = (/o)~^ ° do- 

(3) Compute half the tension of go (i.e., | 5 fQ(a;)|) in each edge of Pi, getting wq e W(Pi). 

(4) Push forward wo to a function vi e >V(P 2 ) by setting, for e e Edge(P 2 ) and y ^ e, 

vi{y) = Y 
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This is independent of the choice of y since g is harmonic. Then retnrn to Step (1), 
using vi instead of Vq. 

Schematically, we are iterating around the following cycle: 


^(r2) 

(2) 

^(ri) 


( 1 ) 


(3) 


w(r2) 

'(4) 

w(ri). 


After each iteration, we can compute the embedding energy: 

Emb(V’i) = rnax 

esEdge(r2) Vi{e) 

Conjecture 6.28. The algorithm above converges to a map with lowest embedding energy. 


Once the combinatorics of the graph map have settled down, this maps in this iteration 
become linear, and the algorithm reduces to finding the largest eigenvector of a linear system 
by iteration. In practice, the algorithm appears to converge rapidly. 

Note the relation to tightness in Equation (1.8): we simultaneously find a representative 
-0 e [0], a map / from r 2 to a metric graph K, and a set of widths tc on Ei, all multiplicative 
on the nose: 


Emb('0) = 


Dir[/ o -0] 
Dir[/] 


EL['0*m;] 

ELH 


where we use a natural extension of EL to a function on >V(r). 


6.7. Extensions and questions. 

Conjecture 6.29. For any strict conformal embedding 0 : Ei ^ E 2 , there is a cover 
0 : El —» S 2 so that [0] is conformally loose. 

Question 6.30. What happens if we vary the definition of SF for surfaces? For instance, we 
restricted to simple curves in Dehnition 6.2. What happens if we drop that restriction, and 
look at general curves? What if we look at the expansion factor of Dirichlet energy for maps 
to graphs instead? What if we look at width graphs, as in Problem 5.20? 

Problem 6.31. Give a direct expression for (some version of) the stretch factor SF[0] or 
lifted stretch factor SF[0] of a conformal surface embedding, analogous to Theorem 3 for 
the case when 0 is not homotopic to a conformal embedding or to Theorem 6 for the case of 
graphs. 

Problem 6.32. In Definition 5.17, we defined a notion of extremal length starting from an 
elastic graph T embedded as a hlling subset of a surface. Give a direct expression for the 
stretch factor (maximal ratio of extremal lengths) between two graphs Fi and F 2 with filling 
embeddings in the same surface E. (Theorem 6 handles the case when F 2 is a spine.) 

Since filling graphs can be used to approximate surfaces, a solution to Problem 6.32 would 
presumably be helpful in answering Problem 6.31. 

Problem 6.33. Extend Problem 6.32 to the more general setting of a group G and surjective 
maps vri(Pi) ^ G and 7ri(P2) ^ G, as in Definition 5.18. 
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7. Dynamics 

7.1. Iterating covers. We finally turn to the dynamical picture: what happens when we 
iterate a map from a conformal surface or elastic graph to itself? Our setting is not quite 
the usual dynamical picture: we are “iterating” virtual endomorphisms, which are not maps 
from a space X to itself, but maps from a cover X to X. 

Definition 7.1. Let Xq and Xi be topological spaces, vr : Xi ^ Xq be a covering map of 
degree d, and (j): Xi ^ Xq be a continuous map. We call this data a virtual endomorphism 
of Xq, also called a topological automaton by Nekrashevych [NekOS, Nekl4]. Dehne to 
be the fc-fold product of Xi with itself over Xq using the two maps tt and 0; e.g., Xq is the 
pullback of the diagram below. 



Concretely, define 

Xk ■= { (xi,..., Xfc) e (Xi)^ I (f){xi) = n{xi+i) } 

Xk comes with two natural maps to Xq: 

• The map tt^ is the map to the leftmost copy of Xq: 

Tikixi, ...,Xk) := n{xi). 

It is a covering map of degree d^. 

• The map (pk is the map to the rightmost factor of Xq: 

It is a composition of lifts of 0 to various covers of Xq. 

We will also use the map 0 to represent the entire virtual endomorphism. 

This construction makes sense even if tt is not a covering map; in this generality, we are 
composing topological correspondences. 

7.2. Asymptotic stretch factors. Now consider the case that X is either a conformal 
surface (S,a;) or an elastic graph (r,a), with a virtual endomorphism 7r,0 : Xi ^ X as 
above. If X is a conformal surface, suppose also that 0 is a topological embedding. Note 
that Xk inherits the structure of a conformal surface or elastic graph from X via the covering 
map TTfc. 

Definition 7.2. The asymptotic stretch factor of the virtual endomorphism is 

(7.1) SFW:-JtoSF[.(.„]'l“ 

Lemma 7.3. The limit in Equation (7.1) exists. 
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Proof sketch. For graphs, sub-multiplicativity of SF (Proposition 6.3) and good behavior 
nnder covers (Proposition 6.22) show that 

SF[cj)k+i] < SF[ 0 fc] • SFM, 

which implies that the sequence SF[0fc]^/^ converges. (It is close to a decreasing sequence.) 
For surfaces, SF does not behave as well under covers, but we still have that if 0 is a cover 
of 0, then SF[0] < SF[0] ^ max(l, SF[0]) by Lemmas 6.9, 6.10, and 6.11; this is enough to 
show convergence. □ 

Lemma 7.4. The asymptotic stretch factor SF is independent of the conformal structure u or 
elastic structure a used to define it. More generally, if f : T ^ T' is a homotopy equivalence 
of elastic graphs with homotopy inverse g -.V then 

(7.2) SF[/o0o^i] = SF[0], 

where 5^1 : F'^ ^ Fi is the lift of g to the cover. 

Proof sketch. Consider Equation (7.2) for graphs and let 0 ' = / o 0 o gi. We have [0(j] = 
[/ o 0n o 5 'n] • Therefore 

SFM < (SF[/]SF[g„])SF[,^„], 

with a similar inequality the other way. Since SF[ 5 f„] = SF[ 5 f] (Proposition 6.22), as n —> 00 
the contribution of the factor SF[/] SF[ 5 f„] to the limit goes to 1, so the asymptotic stretch 
factors are equal. The case of surfaces is similar. □ 

Thus we may speak about the asymptotic stretch factor of a virtual endomorphism 
TT, 0 : Xi ^ X where X is a topological graph or a surface, without reference to the conformal 
or elastic structure. When the covering map is trivial (i.e., the virtual endomorphism is an 
ordinary endomorphism), this recovers W. Thurston’s theory of pseudo-Anosov maps. 

Proposition 7.5. 7/0 : S O A a pseudo-Anosov self-homeomorphism of a surface S (possibly 
with boundary), then SF[0] is the pseudo-Anosov constant of (f), i.e., the exponential of the 
translation distance of the induced map on Teichmuller space. 

Proof sketch. Follows from Theorem 3. □ 

Proposition 7.6. Let T be a ribbon graph, and /ef vr, 0 : Fi ^ F be a virtual endomorphism 
ofV so that 0 extends to a topological embedding X0 : XFi XF. Then 

SF[0] = SF[X0]. 

Proof sketch. Proposition 6.23 relates stretch factors on graphs and on surfaces. The errors 
in the estimates disappear in the limit defining SF, as in Lemma 7.4. □ 

Proposition 7.7. Let T be a graph and 7r,0 : Fi ^ F be a virtual endomorphism ofT. Then 
SF[0] < SF[0]. 

Proof. Immediate from Proposition 6.3, Proposition 6.22, and the definition of SF. □ 

Proposition 7.8. Let F be a conformal surface and tt, 0 : Si ^ S fte a virtual endomorphism 
of F with 0 a conformal annular embedding. Then SF[0] < 1. 
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Proof, version 1. By Theorem 5, SF[(;/)] < 1. Then we have 

SF[0J = SF[0„_i o0] ^ SF[0„_i]SF[0] 

(where cj ): ^ -^n-i is a cover of 0), so 

SF[,^„] < (SF[,#,])" 

SF[^] s; SFM, □ 

Proof, version 2. Here we avoid Theorem 5. 

Let T(0) be the Julia set of 0: the intersection of the images of This set has 

measure 0. Suppose for simplicity that E is connected; an Euler characteristic argument 
shows that E is planar. Pick a conformal embedding of E in CP^. For each x e J(0) and 
n sufficiently large, cfnifPn) can be translated in E c CP^ to a map that misses x. By 
compactness of for sufficiently large n the homotopy class [0„] is conformally loose, so 
SF[0„] < 1 by Proposition 6.14. This then implies SF[0] < 1 as above. □ 

Proposition 7.9. For cj) a virtual endomorphism of either an elastic graph or a conformal 
surface, z/SF[(;/)] < 1, then for all sufficiently large n, SF[(;/)„] < 1. 

Loosely speaking. Proposition 7.8 says that SF[0] detects conformal embeddings in covers. 
Since SF is also independent of the elastic or conformal structure used to dehne the stretch 
factor, this becomes a useful tool for studying rational maps. 

Example 2.7 gives a case where the stretch factor necessarily drops under iteration, so 
SF[0] < SF[0]. 

7.3. Comparison to Lipschitz expansion. There is another natural dynamical invariant 
of a virtual endomorphism of a graph; its (asymptotic) Lipschitz expansion, which we now 
compare with the asymptotic stretch factor. 

Definition 7.10. For F a metric graph and 7r,0 : Fi ^ F a virtual endomorphism of F, the 
Lipschitz energy of 0 was dehned in Equation (4.7). The asymptotic expansion of 0 is 

Lip[0] := jim Lip[0jb''. 

This is independent of the metric graph used to dehne it (up to homotopy equivalence), by 
the argument in Lemma 7.4, . 

A metric train track representative for a free group automorphism : Fn ^ Fn is an 
endomorphism 0 of a graph F inducing 0 on tti so that Lip[0”] = Lip(0)”, i.e., so that 
Lip[0] = Lip[0] = Lip(0). Bestvina and Handel’s theory of train tracks for free group 
automorphisms [BH92,Besll] shows that most free group automorphisms have metric train 
track representatives. (More precisely, any irreducible automorphism has a metric train track 
representative.) 

Given the similarity of the dehnitions, one might suspect these two quantities are related. 
Indeed, they are essentially the same in the case of endomorphisms that are not virtual. 

Proposition 7.11. For cj) : F ^ F an endomorphism of a graph that induces an automorphism 
on TTi, we have 


SF[0] = Lip[0]2. 
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Proof sketch. Suppose for simplicity that F is trivalent and has k edges each of length 1. 
Then, comparing Equations (5.3) and (5.5), we see that, for any curve C, 

\^[Cf ^ EL[C'] ^ l[Cf. 

fC 

We deduce, from Equation (4.9) and the definition of stretch factor, that 

(7.3) -Lip[0]^ ^ SF[0] < fcLip[0]^. 

K 

The error factors in these inequalities go away in the limits defining SF and Lip. □ 

Remark 7.12. Proposition 7.11 may be compared to the fact that the absolute value of the 
largest eigenvalue of a linear operator on a hnite-dimensional space equals the asymptotic 
growth rate of the norm of any vector, independent of the choice of norm. In this analogy, 
EL is like the (square of the) norm while ^ is like the norm. 

For a virtual endomorphism 7r,0 : Ti ^ T, the quantities SF and Lip are in general 
different. By Equation (6.15), we have the rather weak inequality 

(7.4) W[0]<SF[0]. 

See Section 8.6 for a comparison of how these two quantities relate to dynamics on S'^. 

Remark 7.13. The reason it is impossible to hnd train-track representatives in cases like 
Example 2.7 is that the same vertex in T appears multiple times in Ti, and it is impossible 
to arrange for all of these vertices to simultaneously obey the train-track requirements. 

7.4. Extensions and questions. There are many questions raised by this theory. First of 
all, in order to cover the general case of rational maps we need to extend the definitions to 
orbifold fundamental groups. 

Problem 7.14. Extend the dehnition of SF to automorphisms of pointed graphs, or to 
virtual endomorphisms of orbifolds. 

Problem 7.15. More generally, dehne and study SF for a virtual endomorphism of an 
arbitrary group. 

We can also ask about properties of the result. 

Question 7.16. Is SF[0] algebraic? 

Problem 7.17. Give an algorithm for computing SF[0], either approximately or exactly. 

8. Rational maps 

Finally, we come to the original goal of this work, the study of branched self-covers of the 
sphere and when they are equivalent to rational maps. 

Definition 8.1. Fix a hnite set P of points in a sphere S'^. A branched self-cover f : 

P) O is a map so that /(P) c P and so that / is a covering map when restricted to 
S‘^\f~^{P) —> S\P. Two branched self-covers are equivalent if they are related by homotopy 
relative to P and by conjugacy. (Without conjugacy, we could not change the set P.) A 
branched self-cover is rational if it is equivalent to a rational map on CP^, which is necessarily 
post-critically finite. (P contains the post-critical set, but may be larger.) 

We will assume there is a branch point in each cycle of P. For rational maps, this implies 
/ is hyperbolic. 
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8.1. Characterizing rational maps conformally. We first give a characterization of 
rational maps in terms of surfaces. 

Definition 8.2. A surface spine for S‘^\P is a surface E (necessarily planar) together with 
a topological embedding i : T, ^ S\P so that complement of the image consists of one 
punctured disk for each point in P, i.e., so that S\P deformation retracts onto f(E). 

Definition 8.3. Let / : P) O be a branched self-cover and S be a surface equipped 

with an embedding i : T, S'^\P. Then the inverse image /“^(S) is the pull-back in the 
diagram 

f-\^) -^-^E 

i 

S^\f-\P) S\P. 

There are two natural maps: 

• An embedding /-^E) S\f-\P) c S\P. 

• A covering map vr : /“^(E) ^ E. Since it is a covering map, /~^(E) inherits all the 
structure of E. 

If E is a surface spine for S'^\P, then the embedding /“^(E) ^ S'^\P is isotopic to a 
topological embedding 0 : /“^(E) —» E. We therefore have a virtual endomorphism of E and 
can apply the technology of Section 7. 

Theorem 8.4. Let f : {S‘^,P) LD be a branched self-cover with a branch point in each cycle. 
Then f is equivalent to a rational map iff there is a conformal surface spine E for S\P 
so that in the corresponding virtual endomorphism, 0 is homotopic to a strict conformal 
embedding /“^(E) E. 

The above characterization of rational maps appears to have been folklore in the community 
for some time. 

Proof sketch. If / : (CP^, P) O is a rational map, then we may take E to be a neighborhood 
of the Julia set J{f) (hlling in all disks that do not contain critical points). Since the dynamics 
is super-attracting to P on the complement of </(/), the inverse image /“^(E) is a smaller 
neighborhood of J(/) than E, and by taking E to be suitably nice (e.g., using level sets of 
the potential function) we can arrange for /“^(E) c E to be a strict conformal embedding. 

The converse direction is a special case of [CPT, Theorem 5.2] or [Wanl4, Theorem 7.1]. 
The technique, quasi-conformal surgery, goes back to Douady and Hubbard [DH85]. □ 

Remark 8.5. In fact. Theorem 8.4 is true if the assumption that the conformal embedding is 
strict is dropped. With that modihcation, it also extends to more general branched self-covers. 
This strengthening is not needed for the main theorem. 

8.2. Characterizing rational maps nsing graphs. In the setting of Theorem 7, we can 
consider the asymptotic stretch factor of the virtual endomorphism. Theorem 8.4 and 
Propositions 7.8 and 7.9 then tell us the following. 

Proposition 8.6. Let f : {S^, P) LD be a branched self-cover with a branch point in each 
cycle in P and let vr, 0 ; Ei ^ E 6e the corresponding virtual endomorphism. Then f is 
equivalent to a rational map iff SF[(j)] < 1. 
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By the equality of the asymptotic stretch factors for graphs and for surfaces (Proposition 7.6), 
we immediately get the following more precise version of Theorem 1. 

Theorem 8. Let [f] : (S'^,P) O be a branched self-cover of the sphere relative to a finite 
number of points P c S^, with a branch point in each cycle in P. For P any spine of S‘^\P, 
let TTr, 0r : Pi ^ P be the corresponding virtual endomorphism. Then the following conditions 
are equivalent. 

(1) The branched self-cover f is equivalent to a rational map. 

(2) For any spine P for S\P, we have SF[0r] < 1- 

(3) There is some elastic graph spine P for S^P and some integer n > 0 so that 
Emb[0r,n] < 1- 

(4) For every elastic graph spine P for S^P and every sufficiently large n, we have 
Emb[0r,n] < 1- 

Proof sketch. The proof has essentially already been given earlier in the paper. The diagram 
below summarizes the chain of implications between (1), (2), and (3). Here, S is a (conformal) 
surface spine for S‘^\P and 0s is the corresponding surface virtual endomorphism. There are 
various side conditions for which you should see the referenced theorems and propositions. 


/ rational 

A 

Thm. 8.4 


3S: f-\T.) strictly 
conformally embeds in S 


Thm. 4 


3S: SF[0s] < 1 

Prop. 7.8 )£e below 


SF[,#,e] < 1 


3P,n: Emb[0r,n] < 1 

V 

Thm. 6 

' 

3P,n: SF[0r,n] < 1 


Prop. 7.7 


Prop. 7.9 


Prop. 7.6 


^ SF[0r] < 1 


Lemma 7.4 says that the asymptotic stretch factor SF is independent of the conformal 
structure on S or elastic structure on F, giving the equivalence to statement (4). 

The implication 

SF[0s] < 1 ^ 3S : SF[0s] < 1 
is not direct. Proposition 7.9 gives the implication 

SF[0s]<1^3S,n:SF[0s,n]<l, 

which by Theorem 8.4 implies that /°" is equivalent to a rational map. To then conclude 
that / is a rational map and that we don’t need to iterate to hnd a surface spine E with 
SF[0s] < 1, we need to know a little bit more about the geometry. For instance, it suffices 
to know that, except in Lattes examples, some power of the backwards iteration map on 
Teichmiiller space is contracting. This is the only place that this argument relies on the 
original characterization of rational maps [DH93]. □ 
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The virtual endomorphism 0 : /“^(T) —> T is going in essentially the reverse direction to /. 
While / is expanding on the Julia set, 0 is contracting (in an appropriate sense). Therefore, 
for / a branched self-cover with a branch point in each cycle and 0 the associated virtnal 
endomorphism of a spine, we dehne 

(8.1) SF[/]:=SF[0]-0 

This quantity is greater than one when / is rational. 

8.3. Polynomials. One important special case is that of polynomials, where Theorem 8 
recovers known results, which we briefly review. 

Definition 8.7. A topological polynomial is a branched self-cover [/] : (A^, P) O of degree d 
with one hxed point oo e P which is branched of degree d. If snch a map is equivalent to a 
rational map, the rational map is a polynomial. 

Definition 8.8. For a post-critically hnite polynomial p with post-critical set Pu {oo} c CP^, 
the filled Julia set K{p) is the set of points bonnded nnder iteration. It is the union of the 
Julia set J{p) and the attracting basin around each hnite Fatou point, points that eventually 
maps to a periodic cycle with a branch point. The Hubbard tree Tp c K{p) is essentially the 
spanning tree of P within K{p). When Tp intersects a Fatou component, it is required to be 
a union of rays in the canonical coordinates. See [DH85] . 

Now dehne the Hubbard graph Gp c J(p) to be the nnion of 

• the bonndary of the Fatou component containing each Faton point in P, 

• the Jnlia points of P, and 

• the minimal tree connecting the above points in K{p), taking nnions of rays in 
coordinates of Fatou components as before. 

Tp is forward invariant nnder the polynomial p. As such, we can consider dynamics on it. 
To make it a complete invariant of the map, we can decorate Tp with some additional angle 
data at Faton points [DH84, Section 6.1.2]; this is close to what is needed to reconstruct the 
graph Gp from the map on Tp. A abstract Hubbard tree is a tree with an endomorphism, 
marked with this additional data. 

Definition 8.9. An endomorphism / of a Hnbbard tree is Julia expanding if it admits a 
metric for which / does not decrease distances and strictly increases the distance between 
any two Jnlia points. 

The Hubbard graph Gp, on the other hand, is not in general forward invariant, since near a 
point in f~^{P)\P the graph Gp may go through a Fatou component. It does have a natural 
virtual endomorphism, which can be reconstructed from the abstract Hubbard tree. 

Lemma 8.10. An endomorphism of an abstract Hubbard tree is Julia expanding iff the 
corresponding virtual endomorphism of the abstract Hubbard graph has embedding energy less 
than 1. 

Theorem 8.11 (Poirier [Poi09, PoilOj). An abstract Hubbard tree Tp is realizable as the 
Hubbard tree of a polynomial iff it is Julia expanding. 

In the presence of expanding dynamics, we can consider the entropy. 

Definition 8.12. The core entropy h{f) of a post-critically finite polynomial / is the 
topological entropy of the map on the Hubbard tree Tp. 
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Figure 8. Spines for the dendritic polynomial + i 


Proposition 8.13. Let f be a post-eritically finite quadratic polynomial. If f is dendritic 
(i.e., each critical point is strictly pre-periodic), then 

logSF[/] = h{f). 

The statement uses an extension of the theory of stretch factors, etc., to the non-hyperbolic 
case by using marked orbifold points, as mentioned in Problem 8.22 below. 

Example 8.14. Consider the polynomial f{z) = z"^ + i. The critical portrait is 


( 2 ) 

0i—— >i\ -> — 1-I-i^ —z 00 

A spine T and its inverse image T for this polynomial are shown in Figure 8. In the extended 
theory with marked orbifold points, the small loops in Figure 8 should be thought of as 
being of length 0, and ignored in the computation of the stretch factor. If we do this and set 
A sa 1.52 to be the positive root of = 2 -1- A, then the stretch factor of the map 0 : F —> F 
is visibly 1/A. This polynomial equation is exactly the same equation one solves to hnd the 
core entropy. 


If / is hyperbolic (i.e., each critical point ends up in a cycle with at least one critical point) 
and the critical cycles have degree di and length n, (for i = 1,... ,k) then 

logSF[/] ^ max . 

This last inequality is often an equality and often gives answers larger than h{f). 

The fact that / is realized as a polynomial iff the entropy on a topological Hubbard graph 
is positive is essentially a corollary of Poirier’s Theorem 8.11. 


8.4. Comparison to annular obstruction. As mentioned in the introduction, W. Thurs¬ 
ton in 1982 gave a different characterization of rational maps among topological branched 
self-covers, as explained by Douady and Hubbard. Instead of finding an object that guarantees 
that the map is rational, he hnds an obstruction that guarantees it is not rational. We can 
write down that obstruction in our language as follows. 

Definition 8.15. An annular elastic graph is one that consists only of circle components. 
Let A be an annular elastic graph embedded in an orbifold S. Then the join of A is the 
elastic graph Join(A) obtained by 
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• deleting all components of A that bound a disk with at most one orbifold point, and 
then 

• merging all parallel components of A by harmonically adding their elastic constants. 
Harmonically adding the elastic constants means taking two parallel components Oi and 02 
and replacing them with a new component 03 with 

1 _ 1 1 

0(03) a { ai ) 0(02) 

This corresponds to stitching together two side-by-side conformal annuli in the most efficient 
way possible. 

Theorem 8.16 (W. Thurston, Douady-Hubbard [DH93]). Let f : {S'^,P) LL he a topologieal 
branched self-cover that is not a Lattes example. Let Sj- be the orbifold of f. Then f 
is equivalent to a rational map iff there is no annular elastic graph A in Sj and map 
(j): A ^ Join(/“^(y4)) compatible with the maps to S\P so that Emb(0) < 1. 

Remark 8.17. The usual formulation of Theorem 8.16 refers to the maximum eigenvalue of 
a matrix constructed out of A considered as a multi-curve (with no extra structure). The 
above formulation is easily seen to be equivalent. 

Where our Theorem 8 looks for a strictly loosening map 

r'(r)-r, 

the older Theorem 8.16 looks for a (not necessarily strictly) loosening map the other direction 

Join(/“^(7l)) ^ A. 

There is an easy argument that both conditions cannot simultaneously hold (as implied 
by the theorems). If both loosening maps existed, scale A so that Emb [74 ^ T] = 1 , and 
consider the square of loosening maps 


( 8 . 2 ) 


Join(/ 1(71)) ^- A 

rnn ->r. 


The map f~^(r) T is strictly loosening but the map A T is not, a contradiction. 

It is also worth noting that Theorem 8 is technically easier than the existing proof of 
Theorem 8.16. The key analytic point is contained in Theorem 8.4, which is relatively soft. 

8.5. Comparison to domination of weighted arc diagrams. Suppose we are given a 
virtual endomorphism of surfaces tt, 0 : Si ^ Sq, with vr a covering map and (j) a topological 
embedding. Then the associated dynamical Teichmuller space is 

Teich( 7 r, 0) := | S' e Teich(So) | 7 r*S' conformally embeds in S in the homotopy class [0] }. 

Here, Teich(So) is the finite-dimensional space of conformal structures on the interior of Sq. 
(The conformal structures are allowed to have removable singularities or not; that is, the 
corresponding hyperbolic structures are allowed to have parabolic or hyperbolic monodromy 
around the boundary components.) Then Theorem 8.4 and the strengthening in Remark 8.5 
say that the virtual endomorphism comes from a rational map iff Teich( 7 r, 0) is non-empty. 
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Jeremy Kahn has studied the topology of Teich( 7 r, 0) in the context of renormalization 
of polynomials, and has results in terms of weighted arc diagrams [Kah06]. A weighted arc 
diagram is a weighted collection X of arcs with endpoints on boundary components of a 
surface S with boundary. If X is hlling (the complementary components are disks), it can be 
thought of as dual to an elastic graph r(X) embedded in E, with one vertex per component 
of E\X and one edge crossing each arc in X. The elastic length of an edge of r(X) is equal 
to the weight on the arc it crosses. 

Now suppose that you have two surfaces, Ei c Eq. There if X is a weighted arc diagram 
on El and K is a weighted arc diagram on Eq, there is a notion of X dominating Y, written 
X —oY. We do not repeat the dehnition here (see [Kah06, Section 3.6]), but, in the hlling 
case, X do m i n ates Y iff there is a map (j) : r(X) —^ r(K) with Emb(0) < 1 that is compatible 
up to homotopy with the inclusion Ei c Eq. 

Kahn used this notion of domination to control certain inhnitely renormalizable polynomials 
(studying how conformal structures can degenerate). In terms of Teichmiiller spaces, we have 
the following: 

Proposition 8.18 (Kahn, personal communication). //Teich(7r, 0) is not compact, then 
there is a weighted arc diagram X on Eq so that tt*X X. 

Conjecture 8.19 (Kahn, personal communication). If there is a weighted arc diagram X 
on Eq so that 'k*X ^ X, then Teich( 7 r, 0) is not compact. 

We can prove a weaker version of Conjecture 8.19: 

Corollary 8.20. If there is a weighted arc diagram X on Eq so that 7r*X strictly dominates 
X and X fills Eq, then Teich( 7 r, f) is not compact. 

Proof. The conditions are equivalent to saying that r(X) is a spine for Eq and Emb[' 7 r*r(X) ^ 
r(X)] < 1. Then Theorem 7 says that Xir(X) e Teich(7r,0) for t sufficiently small; this 
gives an explicit sequence exhibiting the non-compactness of Teich(7r, 0). □ 

Corollary 8.21. If the virtual endomorphism tt, 0 : Ei —> E comes from a topological branched 
self-cover f : (S'^, P) O with at least one branch point in each cycle, then Teich( 7 r„, 0„) is 
not compact for n sufficiently large. 

Proof. Theorem 8 says that for n large enough, there is an elastic graph T so that Emb[/~"'(r) - 
T] < 1 , which by Corollary 8.20 implies Teich( 7 r„, 0„) is not compact. □ 

In fact, the techniques can be strengthened to show that as long as Teich( 7 r, 0) is more 
than a single point, then Teich( 7 r„, 0„) is eventually non-compact. 

In addition, Kahn has communicated examples where Teich(7r, 0) is compact (but more 
than one point), thus showing that iteration is necessary in Theorem 8. 

8 .6. Comparison to Lipschitz expansion. In Section 7.3 we remarked that there are 
two analytic quantities associated to a virtual endomorphism of a graph: the asymptotic 
stretch factor SF and the asymptotic Lipschitz constant Lip. In the setting of Theorem 8, 
what are the virtual endomorphisms for which Lip[/“^(r) ^ T] < 1? Since Lip[0] ^ SF[0] 
(Equation (7.4)), this is a weaker condition. In fact, Lip[0] < 1 can be taken as a dehnition of 
when associated virtual endomorphism of 'Ki{S‘^\P) is contracting in the sense of Nekrashevych 
[NekOS, Nekll, Nekl4] . 

Bonk and Meyer [BMIO] and Haissinsky and Pilgrim [HP12] studied essentially this 
condition, in the opposite situation to the present paper: they considered the case when 
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there are no branch points in any cycle in P. (For rational maps, this means that the Julia 
set is the whole sphere.) In our language, they showed more-or-less that for a branched 
self-cover / with no branch points in periodic cycles with associated virtual endomorphism 0, 
the stretch factor Lip[0] < 1 iff / has a topological representative as a map on that is 
uniformly expanding. Any rational map has such a uniformly expanding representative, and 
many other maps do as well. 

8.7. Extensions and qnestions. The hrst obvious question is to extend Theorem 8 away 
from the restricted setting. 

Problem 8.22. Extend Theorem 8 to general topological branched self-covers, dropping the 
condition that there be a critical point in each cycle. 

It appears that the theory will extend without much problem to the much larger family of 
maps which have at least one critical point in one cycle in P. (For rational maps, these are 
maps where the Julia set is not the entire sphere.) This requires extending the theory for 
graphs and surfaces to allow dealing with marked orbifold points, as briefly mentioned in 
Section 3.4. 

For the remaining case, when there are no branch points in any cycles in P and the Julia 
set is the whole sphere if the map is rational, there is no clear analogue of the conformal 
characterization in Theorem 8.4 and so it is not clear what the statement would be. 

Question 8.23. In Theorem 8, does n = 1 suffice? That is, for / as given there, is there 
always an elastic graph T so that Emb[/“^(r) ^ T] < 1? 

The answer to Question 8.23 is “yes” for polynomials and for many examples, and The¬ 
orem 8.4 is a corresponding statement for conformal surfaces. Nevertheless, the answer to 
Question 8.23 is “no” in general. See Section 8.5. 

Problem 8.24. In Theorem 8, how does the minimal necessary n grow, as a function of the 
degree and the size of PI 

Problem 8.25. Study the algorithmic complexity of hnding a strictly loosening map 0 : 
/“"■(T) ^ T when / is equivalent to a rational map. In particular, does this give a polynomial¬ 
sized certihcate that / is rational? 

Presumably a solution to Problem 8.25 would hrst require a solution to Problem 8.24. 

Problem 8.26. Study SF[/] as / varies over rational maps. How does it compare to known 
invariants? 

Question 7.16 is one starting question. The connection between SF and the entropy on 
the Hubbard graph described in Section 8.3 is relevant, but in general there is no reasonable 
forward dynamics on P. 

Finally, Theorems 8 and 8.16 give opposite combinatorial conditions. These should 
presumably be related. 

Problem 8.27. Prove combinatorially that, for any virtual endomorphism 0 : P ^ F of a 
graph, either SF[0] < 1, or there is an annular obstruction in the sense of Theorem 8.16. 

A solution to Problem 8.27 would give an alternate proof of the main Theorem 8, using 
W. Thurston’s Theorem 8.16. 


FROM RUBBER BANDS TO RATIONAL MAPS 


47 


Another family of problems comes from applying Theorem 8: Prove that certain maps 
are rational by finding explicit maps 0 : /“^(P) ^ P with Emb(0) < 1. We give one sample 
problem. 

Problem 8.28. Rees and Tan Lei [Tan92], and later Shishiknra [ShiOO] showed that two 
quadratic polynomials are topologically mateable iff they do not lie in conjugate limbs of the 
Mandelbrot set. Reprove this with the techniques of this paper. (See Examples 2.4 and 2.5 
for one approach.) More generally, find a criterion for cubic or higher-degree polynomials to 
be mateable. This is likely to be more difficult, given the cubic example of Shishiknra and 
Tan Lei without Levy cycles [STOO]. 

Warning 8.29. A topological mating is not the naive (formal) mating obtained by gluing two 
polynomials at infinity. See the survey by Buff et al. [BEK+12]. 
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